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ABSTRACT 


The 9catteting of phonons from a substittitioHai inpT*“rity in 
cubic lattices is studied in detail with the help of Gieen function 
technique used by Lifshitz and Montroll end Potts. Exact solutions 
for scattering amplitude are found explicitly and the scattering 
amplitude is expanded in terms of partial waves characterised by the 
irreducible representations of the point group of the lattice. The 
scattering cross sections are found in the limit of long waves showing 
Rayleigh type of scattering in absence of resonances which are taken 
care of 1:^ the reaoriance denominators. 

Specific exanples of the acatteiing of phonons from suhstltujfcional 
iiipurities in simple cubic, body centred cubic and face centred cubic 
lattices are worlced out in detail* Cunres are given showing the 
dependeme of resonance frequencies on parameters characterizing 
changes in mass and force constant* The long wavelength behavior of 
the scattering process is discussed and the total cross sections are 
found using optical .‘theorem in this limit • The dependence of the 
scattering cross sections on the direction of incidence is shown by 
noting the different partial wave contributions for different directions 
of incidence* The inverse of mean free path for phonons due to 
scattering by substitutional imparities is found in the long wavelength 
limit# It xs shown to be proportional toi^ ^ with resonance 
possibilities and the interference terms between partial waves 
characterized by different irxeduciblje representationsare given 
eDsplicitly# Phonon scattering from a substitional inpurity in a 



diatomic einple cubic lattice is also discussed and clirves are draim 
showing the dependence of resonance frequencies upon the parameters 
characterizing changes in mass and force constant* 

At the end, a model of body centred cubic lattice is consideied 
in lAiich the polarization of phonons is taken into account in the 
decciiption of the scattering process, which is shown to depend upon 
the direction of incidence as well as the state of" polarization of the 
incident phonons. The dispersion law for the model does not go over 
to the Debye model in the limit of long waves and consequently the 
dependence of the number of partial waves contributing to the scattering 
cross section even in the limit of long waves depends upon the incident 
direction and polarization of incident phonons. The result for the 
isotope defect case in the limit of long xraves differs from that 
calculated by assuming Debye model by a numerical factor which depends 
upon the direction of incidence* A compariso,^ is made between tbis 
and the so called scalar models, • ; 
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CHAPTER I 


IFTROITJCTIC^J 

The problem of scattering of phonons by e substitutional 

1^2 3 

impurity was first treated by I.M. Lifshitz whose pioneering work 

on the problem of Dyhhmic Theory of Won^ideal crystal lattice is well 

4>-»1 *7 

known. It has subsequently been discussed by maiy workers including 
-| Q— 19 

the author . Most of the fundamental concepts and methodological 

aspects of the problem are present in Lifshitz’ s work. . It has been 
shown that the problem can be worked out exactly using Gre«n function 
techniques. The perturbation due to impurity may be large but the range 
of this perturbation is assumed to be quite small within the crystal# The 
principal result is that the perturbed region surrounding the impurity 
appears as a separate system whose energy levels depend upon the energy 
spectrum of the host lattice. If the perturbation is sufficiently large, 
an inprrity level may lie outside the band, xd.th a long lifetime since 
there are no states in the neighbourhood into which it can decay. On the 
other hand for moderate perturbations the impurity level vdll lie inside 
the band with a finite probability of decay. The latter leyel has 
profound effects on the scattering of phonons from the iitfjurity. • 

Is an exHuple Lifi^itz, considered the case of an isotope defect 
and the scattering of phonons from it. He obtained the expression for the 
acuttecring ampli'ttde by expanding the proper Green function for large 
distances. The behavior of this process in the limit of long waves arc 
also discussed. Montroll and Potts independently 'Considered the 
problem in the same way as Lifshitz did. They also showed that the 
• case of diatomic lattices can be handled by the formalism for mono- 

atomic lattices by performing the so called M* transformation .Although 
Lifshitz noticed the fact that by the infcroducticsn of a substitutional 

« *1 • • 



impurity the translational symmetry of the lattice is destroyed and 
only the point group symmetry remains, the idea of expanding the 
phonon scattering amplitude in terms of partial waves characterised 
ty the irreducible representations of the point group of the lattice 
was developed by others 

The rest of the work in this field has been mostly matters of 
detail. The literature is extensive hit the notations are uncorrelated, 
■some times to the extent of confusing the reader. Maradudin*^ has 
discussed the details and refinements of the problem on several 
occasions. Takenq^'^, Callaway®, Wagner^^,Kiumhansl^^'”^’^ and Thoma and 
Ludwig have also considered the details of the problem with its 
implications on the e'valuation of thermal conductivity due to scattering 
of phonons from substi-fcutional impurit^.es , In many of these the ides 
of resonant scattering has been stressed. Krumhansl and Matthew^*^ 
considered the problem in a one dimensional model in a rather complete 
way, and have demonstrated interference between different partial 
waves in the expression for relaxation time. This has been elaborated 
for three dimensions in the present work. Elliot and Taylor^^ solved 
the problem by using double time Green function technique in an 
elegant way. Recently Klein has given a completely general 

formulation of the problem within the franewoit of the quantum 
mechanical formal theory of scattering. Several important aspects 
like thfe optical theorem end orthonormality of the solutions have been 
discussed in this work, Wagner^^ has discussed the case of molecular 
inpurlties by introducing the molecular Green function. 

There has been considerable overlaping among the appieaohes 



and results of the vaiious workers mentioned above. The overell 
picture at present is that we have made substantial progress in 
understanding the details of the problem, but this by no means is 
complete and final. Experiments in this field has lagged behind theory 
and new experiments may reveal more details than has been known 
hitherto. The main source of experimental information on scattering 
of phonons from substitutional inpurities at present is the low 
tenpereture thermal conductivity data. By analysing this. Walker and 
Pohl^^“^2 were able to infer the presence of resonance scattering 
leading to the deviation from KLemens' result^^. Baumann^^ and KLein^^ 
also arrived at similar conclusions . Very recently Huebener^'^ has 
estimated the scattering cross section for phonons scattered from 
vacancies by measuring the change in thermoelectric power in Platinum. 
A remark relevant to these experiments is that usually the inpurities 
do not always go into substitutional sites. They sometimes go into 
interstitial positions so that the results have to be interpreted 
cautiously. 

There is a great deal of similarity among the problems 
• of scattering of phonons, band electrons and spin wraves from 
substitutional impurities. This has been explored by Callaway® in 
detail. Formalisms identical to that for the phonon scattering for 
the other two cases have been developed by various works . 

In this dissertation scattering of phonons from substituti— > 
onal ImpuTities in some three dimensional lattice models are conside- 
red in detail. The scattering airplitudes are expanded in terms of 
partial wayes chaiacterised hy different irreducible representations 
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the point group of the lattice anri the scattering cross sections 
are evaluatec^ in the limit of long waves. The question of resonant 
scattering of phonons in the optical banr^ of a diatonic simple cubic 
lattice is discussed. At the end a model of body centred cubic lattice 
IS considered in which the polarisation of phonons is taken into 
account in the description of the scattering process. The question of 
ths dependence of the scattering process on the direction of incident 
phonons is also examined. Finally the implications of a departure from 
Debye model in the limit of long waves is discussed. 



CHAPTER n 


GEI^trml FORMJUTION 

S 2,1 Introduction . 

The theoiy of scattering of phonons from substitutional 

imouiities started with efforts concentrated on isotopic impurities 

irrrolving only a change’ in mass. It was further believed that the first 

Born approximation results were adequate to d.escribe the process , 

but thio assumption was soon abandoned . The low temperature 

thermal conductivity experiments by Klein^®, Walker and Pohl^^“^^ and 
25 

Bau'cuonn ' showed the inadequacy of first Bom approximation results 
and the necessity of including resonant scattering which oan be done 
by solving the problem exactly. They farther showed that the rather . - 
simple case of isotope defect was also inadequate. One should consider 
the scattering of phonons from point defects characterized by changes 
in maes as well as force constant. Moreover one ^ould look for exact 
solutions rather than perturbation approximations. We consider such 
an aporoech to the problem of scattering of phonons from a suhstitutional 
impurity in the present work. 

Some other investigators have considered this approach as 
already mentioned in the last, chapter. Many other workers^'^^^ have 
been interested in the determination of only the localized modes due 
to a 'general substitutional impurity' by which we will mean, 
throughout this work, a substitutional impurity characterized by 
changes in mass as well as force constants. However this dete imination 
of localized modes can be achieved in a natural way from the solution 
of the general scattering problom. 


« • 
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The simplest mcyiel for the scattering nf phonons from a 
gur.eral substitutional impurity has been consic^erec^ by Kiumhansl and 
Matthew^^. Although they consider a one-dimensional problem, the 
qualitative conclusions derived are quite instructive and useful. 
ICLein^"^ develops a t-matrix formalism which is applicable to three 
dimensional crystals. This has been extended further b/ his latest 
work^^ which contains many important theorems and an averaging 
procedure for a distribution of defects. It is a self consistent 
formalism applicable to many models and yields results that 
identical to those of the simple formalism described below ( when 
applied to the same models), except that it is perhaps easier to 
exLract numbers out of the latter than from the former formalism. 

Also the explicit detailed treatment of the resonance denominators 
in the latter formalism facilitates the investigation of resonances 
and localized modes. 

Before discussing the present formalism it is worthijhile 
to note that this scattering problem is analogous to the potential 
scattering in ordinary quantum mechanics except for two important 
differences. It is a particular case of the more general problem of 
scattering of quasi-particles in crystalline solids 'vdiich has been 
s'cudied by Callaway^. One of the differences is due to the fact that 
the dispersion law giving the wave number dependence of the single 
particle energy reveals the typical band structure involving both 
an upper and a lower bound to the energy. This permits the possibility 
of occurrence of bound states both above and below the optical bands 
and above the acotistic hand, depending on the nature of the 



.. 7 .. 

scatterer. Moreover for the scattering process both the incident and 
the scattered particles have their admissible energies limited to the 
band. The wave number is also limited to lie within the first Brillfxxin 
Zone. The other difference is relate'^ to the symmetry of the problem. 
The interaction due t-'^ the scatterer usually possesses the symmetry 
of the point group of the crystal rather than the spherical ^rmmetiy of 
the potentials studied in the theory of potential scattering. This 
feature necessitates a change in the method of partial wave analysis of 
the scattering amplitude which now has to be decompose'^ in terms of 
Partial characterized by the irreducible representations of the 

point group of the lattice. 

The present formalism is perhaps mathematically simpler and 
less formal than those cited above . The scattering amplitudes and 
their partial wave expansions are explicitly presented and tho optical 
theorem in the usual form is employed to demonstrate the non— central 
effects of a genera^l perturbation. The solution is exact and includes 
resonant scatterings. The long wavelength limit is found easily. The 
interference between terms corresTJonding to different partial waves 
can also be seen explicitly in the evaluation of relaxation times. 

§ The Model 

The formalism developed here depends upon a number of 
assumptions. The first is the adiabatic approximation in which it is 
assumed that the electrons instantenously take up a configuration 
appropriate to that of the displaced nuclei and that the change of 
energy of the electrons in the distorted lattice contributes to an 
effective intemuclear potential. This is quite plausible for the 
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vibrations in the acoustic moies that we will deal with. The other 
important approximation is the harmonic approximation. Since the 
harmonic force constants assumer^ here may have a temperature dependence, 
we are in effect assuming a pseudo-harmonic approximation which is 
adecpJiate for low temperature properties. 

Here we consider a monoatomic lattice and investigate cases 
where the question of polarization does not enter into the picture. The 
so called * scalar models' belong to this categoiy. In addition there 
are models other than these scalar models which do not mix polarization 
such as the one discussed in the next chapter. The discussion of 
polarization is taken up in chapter V of the piesent work and a 
discussion of a diatomic lattice in the next chapter shows how the 
results for the monoatomic lattices can be easily carried over to the 
case of diatomic lattices by using the so called M transformation^. 

With slight modifications; the results obtained here can -be carried 
over to the case of the scattering of spin vmives arid 'the 

scattering of band electimns^*^^^ from substitutional impurities in 
the lattice. 

The lattice under consideration has N atoms of mass M each 
which ere coupled to their neighbours thrmugh harmonic force constants. 
The time independent equations of motion can be written as, 

AU - ( M*/) U = 0 ( 2.1 ) 

A 

where A is the force constant matrix of dimensions ¥ x F 
and U is the N-dlmensional vector whose elements are the one 
conponent displacements ( scalar field ) of atoms from their 
equilibrium sites v&th a time dependence of 1116 form e3cp('-d.a>t ). 
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The matrix A can be shovm to be cyclic and hence its 
eigenvectors Fq elements of the form 

(k )= exp (Ik .R ) 


( 2.2 ) 


whore R denotes the lattice sites with reference to a given 
origin. The vector k ( with R different values) specifies the 
eigenvectors and iiie corresponding eigenvalues. Ro (k ^ S ) kas the 
f<~'Ta of a running wave or a phonon propagating with a wave vector 

k * Substituting e^^u&fcioh (2-2) in equation (2.1), one gets the 

. A 

eigenvalues of A which give the frequency squared as a function of 
k i,e the dispersion law which is written as a) (k ) and this has 
both an upper and a lower bormd defining the band. The distinct 
values of k lie within the firet Brillouin Zone. The eigenvectors 
So ( is ) satisfy the orthogonality and closure properties in an 
obvious way 


and 


n C(k , £0 Uo (k , R ) - g^R 


S' 


23 D„*(k<, S ) Ck , B ) = ?k, k' 

R 


(2.3) 


These properties of A and its eigenvectors are direct 
consequences of the cyclic boundaiy conditions imposed on the 
lattice , 

_...g.3.The perturbed lattice . 

Let a substitutional impurity be put into the lattice . This 
changes the mass of the atom at the substituted site to M + M 
and the obupling of this impurity atom to its nearest neighbours 
is given ly the changed harmonic force constant liiere 

^Ydenotes the change in the normal value Y* . Thus the impurity 
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affects the site it oooapios plus its neighbours uhcso number is 
fleteimined hy the range of interaction, ill these sites including 
that of the impurity atom are osUed the sites affected by the 
perturbation. In most of the cases these affected sites can be assuMd 
to hare the symmetry of the point group of the lattice with the impurUy 
sire as the point left inraidant under the opoiatins of this point 
group. Strictly spcaUng, the resulting symmetry group of the system 
uiU bo the Intersection of the point group „f the lattice and the 
symmetry group of the perturbation. Wha.t w, are assuming here in 
effect amounts to helievlng that the symmetry group of the 

porturbaticn has at least all the elements of the point group of 
tho lattice. 

Then the displacements of all the atoms at the affected sites 
form a representation space of the point group of the lattice in rtilch 
ttc reducible representation of tho point group can be decomposed 
in terms of the irreducible representations p^il) ot the point 
gmup ea 

"here are sere or integers. The aboTe representation space 

wixl be called the ^space ♦ 

The equations of motion of the perturbed lattice can be 
written as 

2 - ( M co^) 2 = P 2 ( 2,5 ) 


tho N X IT pertuihation matrix P has non>. ranishing elements 
only in ^ P*, space and can be partitioned into the form 

I 0 0 


( 2.6 ) 
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whaiG the matrix p is of dimensions n x a and n is the number of 
lattice sites affected by the perturbation, i.e. the din^nsionality of 
the p -space . 

The solution of equation (2,5) outside the unperturbed band 
will give the discrete frequencies corresponding to bound states in 
potential scattering. They are the localised modes of the lattice 
with impurity. The scattering state solutions are obtained for values' 
of u) within the unperturbed band and can be written as satisfying 
an equation of the form , 

( 2 . 7 ) 

where the Green function matrix G is defined as , 

A .A P A 

G=-.(A-McGI) (2.8) 

where ? is the unit matrix, Uq is the unperturbed solution having 
the same eigenvalue as U this is always possible in a large 

lattice^ because the frequencies form a continuum • 

The matrix elements of G can be written explicitly as^ 


G ( J , E1)= _ i V n<,(k ,£ ) u‘(£,£') 

^ <0^(k)- J 

The summation over k is over the first Brillouin Zone. 

The disadvantage of the use of such Green functions for the 
three- dimensional problems is that they must be calculated 


(2.9 ) 


numerically even for the simplest cases if quantitative results are 
desired. However some numerical tables are available^^^^^^^^ end 


recently a simplified numerical method based on the Fourier series 

34: 

expansion has been given by Mahanty , These Green functions also 
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occur in problems of scattering of band electrons ant^ spin waves 
from substitutional impurities. 

Using equation (2.2) for / E ) in .equation (2.9 

the Green function matrix elements can be written as, 


G(^^»S ,£1 ) 



e3cp {_-± £•(£’-£, ) ^ 

cj? (k )- cu^ 


( 2.10 ) 


Ihe matrix G is cyclic like the matrix A and. hence its elements 
depend upon the difference H - g' . 

Now the equation (2.7 ) can be solved for U . The first 
Bom approximation solution consists of writing Uq for U on the 
right hand side of equation (2-7 ), Thus the first Bom approximation 
solution of equation (2.7 ) reads , 

a =210 - S P Eo (2.11 ) 


This shows that U Is equal to the incident wave plus 

another part which we wish to represent an outgoing wave for our 
scattering problem. This can be achieved by writing equation (2.10 ) 
in the form , 


^ ( k) } S’ 


1 

inr 



^ (k ) - 


1 


( 2.12 ) 


with £ -5.0 . 

In the usual way the second Bom aporoximatinn solution 

will be , 


S = 26 


(2.13) 



t 


■.13 


ar/ so on. Thus one gets the perturbation series solution. But this ' 
^oes not show resonance in the scattering process. Also for lai^e 
prrtuibations such a series has convergence problems. 

Instead of equation (2.13 ) we wish to solve equation (2.7 ) 
exactly. In a formal way, the exact solution reads. 


A A , A A A. 


TT T-r I' , /\ O 

2 = 2o - G P( I r G P)- 


U. 


(2.14 ) 


The orthonormality of solutions of the foim oi »<luHlon(2 .14) has 
been discussed in detail by IQein and others'^. We now proceed 
to describe the operational procedure for obtaining an exact 
solution of equation (2.7) like the one given by equation (2,14), 

^ ^ Solution for the scatt e ring amplitiida 

In order to obtain an explicit solution we partition TJ in 
the form, 

- ' (2.15) 


AdJ 


whore u is a n-dimensional vector in the P-space corresponding 
to the nonvanishing partition p of f described in equation 

(2.6). In a similar fashion, the matrix G and the vector W can 
be Partitioned into , 


and 


A 

G 


So = 



( 2.16) 


( 2.17) 


where g is n X n matrix in P^epeee correepondii^ to the 
partition & and 310 is in conaspoMing to a . 
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Let us define the N *- dimensional vectors and Sq 'bgr, 

(2.18) 


S = P TJ 


and 


So =PlIo 


(2.19) 


By multiplying the partitioned fom of P given by equation (2.6 ) 
and using equations (2.15) and (2.17) for the paitition of U and 
Up one gets. 


I ^ 

V2 


and 


So 


£0 

a 


(2.20 ) 
(2.21 ) 


where la = p u and s 


so 


^ He 


(2.22) 


Now multiplying equation (2.7 ) by P on both the sides 
by the left, one gets 


or 


A A /\ A A 

Pu “ PUq -pgpu 


A A 

s » ^ - P G S 

— —o — 


i 2.23 ) 


which in the partitioned form reads. 




0 

0 


whence the reduced equation for is , 

A A 

a “ £0 - p g s 


so that |[ «= (l + pg )*~^ s 


So 


or 


- i» / T ^ ^ V*"l A 

S“(I+pg) pUo 


(21) ^(22) 


G(12) 

A 

G/ 


(2.24 ) 

( 2.25 ) 
( 2,26 ) 
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Thus the matrix irrvolvef? in equation (2.26) is only of the 

. . f\ A, 

iumenflions n x n. For a given case one knows ? and g explicitly. 

But still to fire^ the inverse c'f(l + ?g)ina stiaight-forward 

manner is often veiy tedious. To ease this task, one can employ a 

A A A 

unitary matrix Y vhich block diagonoli 2 es 1 + P g. One could eYStt 

think in terms of a imitory matrix which completely diagonolizes it, 

but such a matrix is difficult to obtain. On the other hand there 

A 

exists a unitary matrix V T»Aiich can be found easily by group 
theoretical procedure and which can block diagonolize p and g and 
hence ( I + p g ) into blocks of quite manageable dimenBions, Then 
it is very easy to find the inverse of( I+pg). The number of 

such blocks, the dimensionality of each block and the number of 

A 

identical blocks for the block diagonolization by such a T matrix can 
be knoi«i by applying the following theorem from group theory to the 
decomposition of P given by equation (2.4 ) . The thoorem^^ states 

that if r and if l~^ has the dimensioneli'ty 

V 

rijf then in the block diagonal form there will be identical 

X a^ blocks. 


Further if one writes the block diagonal matrix 



these blocks which corresponds to a certain irreducible representa- 
tion is non-zero, then using each of them one can evaluate 
vectors s^ whose sum will be s . This deconposition of s ^ ^ 
is equivalent to projecting the different symmetry adapted functions 
out of s , which belcxig to different irreducible representition* of . 
the point group. This in turn leads to the partial wave expansion 
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of the scattering amplitude. Thus the matrix 7 found by g^up 
theoretic procedure not only facilitates the evaluation of the 
inverse matrix in equation (2.26 ) but also leads to partial wave 
expansion of the scattering amplitude as shown below . The evaluation 
of such matrices is described in the Appendix HI. 

Thus let, 



where T is the blor.k diagonal matrix with k blocks. It is veiy 
easy to find the inverse of T since this has KLocko of veiy small 
dimensions. Let T be written in the block form. 



( 2.28 ) 


Now each of the blocks t^ , t 2 , ••• tj^ correspond to a unique 
irreducible representation of the point group of the lattice. ¥e 

define the matrices such that has non-vanishing elements 

in the l) th block only, i.e 



(2.29 ) 


and 



h 


z: 



(2.30 ) 


mi— ^ ^ A A . 

Hius P g ) »= V V 


h 

S 


A A A. 

V r 


^ Z SL 
>‘*=1 
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, ^ A ^ 

where >= V V p ^ 


( 2.31 ) 


One obtains the same result if one block '^iag'^nolizes p as 
A 

well by the V matrix and writes it as the sum of matrices like 

/N. 

f mil ti plies the th blocks of these and defines s,i as 
s jj, = V (m p Ho • 

Th® s_ defined like this belong t'^ particular nrows of 
particular irreducible representations of the point greup of the 
lattice-. This identification can be done by iriGntif 3 d.ng the columns 

A 

of V - jaatrix which is done during its evaluation. One would get 
the same result if one finds s hy direct inversion in equation 
(2.26) and then applies the group theoretic projection operators 
to 3 to find the s ^ that belong to given rows of given 
irreducible representations. 

Now going back to equation (2.20) one finds that S can be 
written as , 


S <= k 

<_.7 

:)>=i 




where S^, 




0 


( 2.32) 




and then equation ( 2.7 ) reads 

h 


>>=1 


A 

G 




( 2.33 ) 


This equation is the analogue of the familiar partial wave 
analysis result of the theory of pot2ntial scattering. Here the 
number of terms is finite and equals the number of submatrices in 
the block diagonal form of p uhich in turn depends upon the range 
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.of inte lection dtie to the impurity end the point group of the lattice. 

In component fom, equation (2.33) becomes, 

h 

U(k, R) =Uq(]c,R) G( S^Sn) (l^R ) (2.34) 

vdiere ^ refers to the lattice sites affected by the perturbation including 
the impuritj" site and the summation is over all of them. These are 
quite small in number and the non-vanishing elements of refer to these 

only. Using equation (2.10) for G( one gets for large R, as 

shown in the Appendix U, 

- - Z1 (fio-(S-40 /I si 

'R-^'bO y. 

—O 

and consequently equation (2.34) gives 

U(^R) = U^(^R) +^2^g (1^) exg ji ko-R)'^ (l^k^) (2.35) 

all ^ fRiir^ “ 


fsp 


I) 


where Fo (k,k ) is a function of the components of k and k . Here kQ 

^ 2 ^ I 

is such that (fej) = cD with (k) I ^ having the direction 

of R. For symmetiy directions in the lattice ^ represents the 

direction of the scattered wave and often there are not mole than one' 
2 !>■ 

satisiying Cd ^i£o) = along these directions. 

Equation (2.35) gives us the desired solution with the outgoing 

spherical iraves characterized by e3cp(i kQ»R)/|R J idiose co-efficient 

is "the scattering anplitude. This point has been discussed in detail 
9 

by Callaway < 


Thus the scattering amplitude is 
h 

f P., (k,k ) 

y=i ^ 


(2.36) 
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an(? the desired partial wave expansion can be written as^ 

where f^(k,kQ) = si'^) with >' summed over values 

where is the dimensionality (degeneracy) of the yi^th irreducible 

representation. 

For incident waves along the symmetiy axes in the lattice, it is 
quite easy to evaluate the forward scattering amplitudes. Thus one sets 
k “ i^Q In equation (2.36) to get the forward scattering amplitude. 

In describing the details of the partial wave analysis we shall 
follow the classification of S,P,D,F etc. "type of waves according to the 
Wigner's scheme described ly Callaway^ and given in detail in Appendix III, 
Thus having obtained the scatteiing amplitude one can evaluate the 
differential scattering cross section in the usioal way. Also the total 
scattering cross section can be found by integrating over the angles. 

Such a procedure is cumbersome in general, but in the limit of long 
waies it can be done for sinrole cases. In this limit the constant 
frequency surfaces become spherical for the model under discussion end 

hence CJe) i lies along k^. Then represents the direction 

I '«*o 

of propagation as vrell as the direction of observation for all directions 
of incident waves which are characterized by k and the angle between k 
and ^ is the scattering angle. 

As mentioned above, one can evaluate the forward scattering 
amplitude for incidence along symmetiy directions . For siich directions 
of incidence one can evaluate the total scattering cross section 
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by usin^' the optical theorem. That the optical theorem holds good for 

thci type of scattering under consideration has been discussed by Thoma 
13 15 

and Ludwig and Klein . ¥e feel that the optical theorem given ty 
13 

Thoma and Ludwig vfcich can be derived using the results for Q — i 6 ) 
given Callawsy , is sufficient for the present work and it is given 


by 


cr: 


^/k)j 

I STBd ^ ^^(k) j 


i 


Im ^ f (k,k) exp^ ( e^-S)^ 

(2.38) 

where f(k,k) is the forward scattering amplitude. Here (k) 

means the second derivatives of cd^(k) Js space with respect to a ' 
coordinate system that is perpendicular to grpd^cd^(k) and diagonalized 
the tensor cO j ..(li)> (ij 0 and is given ty sgn 0)^L,(k). 

However, it must be remarked here that from the experimental point of 
view one is interested not in the total scattering cross section but'the 
momentum transfer cross section which is required for the calculation of 
thermal conductivity. This is discussed in the following section. The 
linri-t of long waves is veiy relevant for such calculations because it 


is a plausible approximation for explaining the results for low temperature 
experiments. Also the optical theorem in this limit takes the usual 

k)j|for the model under discussion, 

Ili is clear that resonant scatterings are included in the above 
calculations. Equation (2.26) shows that s has det | I + p g| in the 
denominator. Further from equation (2.27) it follows that det | I + p ^ | 

= det j T~j which in turn is the product of the deteradnants of the blocks 


form of C57 


47r 

“TT 


Im 


[5(k, 


of ’7“ which correspond to different symmetiy types. Thus det | I + p g| 

- dot I irl = 7r. where is the determinant of the block which 
coiresponds to a given symmetiy. These are exhibited explicitly in the 
nexu twn chapters. These sub-determinants b. occur as denominators in the 

J 

expressions for f^^(k,k^) or equivalently one can speak of TT, b. 

occ-'ring in the denomina^tor of f(k,k^). In general b. are cLplex within 

the band and can be written as b = b + i b. , wherc b. and b 
are real. Moreover these are functions of tho incident eneig;^ E. Outlide 
the band and b^, are completely real. Localized modes appear when 

anj- of the sub-determinants b^ -vanish outside the band . This can be seen by 
multiplying equation (2.5) by $ from the left to "get U = -G ^U. Here 
we do nob have a because we consider outside tho band. Mow follovring 
the partitioning described by equations (2.15), (2.16) and (2.6), we get 

A /s 

- “ “ S P u (2.39) 

This is an eigenvalue, eigenvector problem in p" -space. Solution of this 
equation exists for Ct)^ determined by the condition det | I + g p | = TT b 

■= 0. These are the localized modes of -vibration whose amplitude falls off 
rapidly with distance from the iupurity and these have been discussed in 
detail by many workers^"^^. Since b are identified into different ^mmetiy 
types, one speaks of localized modes of different symmetries and degeneracies 
corresponding to the b^ that vanishes to give rise to these modes. Since 
we evaluate the factors b^. explicitly in the present foimalism, it becomes 
very ea^ to investigate the existence of localized modes. 

Inside the band, the real part b 

hr 


may vanish for certain values 
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of incident eneigy E. It will correspond to a resonance in the scattering 

amplitude provided that the imaginaiy part b. . and ^ '^j^r have the 

* ^ ^ ^ 
same sign and their ratio is small. This additional requirement comes from 

the consideration of the density of phonon states of the perturbed lattice 

as shown below. We do not use the density of phonon states elsewhere 

although it can be evaluated easily using the above formalism. The 

density of phonon states for the perturbed lattice can be written as 
1 . A . 

~ ~ "rr ^“(Tr G’ ) where the perturbed Gieen function G* is given ly 


A A A A A A A _1 

G' = G - G P (I + G P) 


whence the contribution of the impurity to the density of states can be 
written as 


An = - la (Set |l + G p| (f..40) 

A 

But due to the partition of P into the only non-vanishing small block 
P; one has 




1 ^ A A. 1 

j( I + G p : 

)| = det ^ 

I + g P 


j 3 

»here »B also use the fact that (l + p g _ (f + J ) and that the 

transposition does not alter the value of the determinant.. Then equation 
(2e4o) becomes ,, 

dE 






db 


3 ,: 


db. . 

1 >3-: 


dE 


2 ,t dE 


L 




( 2 . 41 ) 


..23 


But at resonance, b =0 and ^n.=. _ "I (2,42) 

3.1* ^ ttL 

is the contribution from a particular impurity state, which when 
compared with the resonance fomula 

1 Pw 


n 


(E) 


TT (e-E )^ + p ^ 

'• O'' lu) 


which at resonance reads 


n(E) = 


TT fl 


Hi 


gives the result for the level width 


11 

CO 


- 


bl 


(2.43) 


Here prime denotes the derivative with respect to E. For resonance 
Ho should be positive. The sm8.11er the the sharper is the 

rosonanco . 


^2.5 Relation with the transport phenomena: 


So far only one defect at a particular lattice site has been 
considered. In actual practice one has a number of such defects randomly 
distributed over the lattice sites. These can even migrate and interpret 
with each other. But at low temperatures and small concentrations, it is 
perhaps not a bad assumption that veiy little migration of defects takes 
place. We want to consider the effect of scattering under these conditions 
upon the heat transport phenomenan^ but since we will not be interested 
in the explicit evaluation of thermal conductivity in the present work, 
the details will be omitted . 

The problem of thermal conductivity can be tackled in two ways. 
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The CO iventional method is that of wiiting down the Bolt 2 imann transport 
XJ 36,37 

equation and solving it in the Kinetic theoiy approximation* The 

other more recent method is through the so caU^d Kubo^® foiraula. This 

39 

has been used recently ty Behera and Deo among others to calculate the 
tneimal conductivity. However we confine ourselves to a discussion of 
the cox"vent5-onal method. 

In the conventional method one has two distinct approaches which 
give the same result. One is the classical method in which one warites 
tue energy flow eejuation along the lines indicated by Choquard^^, evaluating 
the various quantities for the perturbed ^stem. The other is the 
quantum approach through transition probabilities and scattering aross 
sections. Since in the present formalism one can evaluate the scattering 
cross sections, we discuss here this latter approach. One starts with a 
Boxtzmann equation for the n’umber of phonons with momentum k and 
studies its time evolution through the equation of the form 



During thermal conduction, the left hand side is due to the temperature 
gradient and has the form 

^.grad T (-^) 

whereas the right hand side is due to the scattering for the present 
case. In the relaxation time approximation one writes 



(2 .45) 
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which shows that in absence of the temperature gradient, the deviation 
fi-om equilibrium in the mode k damps cut esiponentially in the charaet eristic 
tine (k) . With this one can solve the Boltzmann equation (2.44) for 
and cf’lculate the thermal conductivity. In terms of heat capacity 
G(k) and the group velocity y, the expression for thermal conductivity is 


K 


V 


(2Tr) 


J d\ t;(k) C(k) v^ cos^0 


(2.46) 


where 7 is the volume, 0 is the angle between y and grad T and the 
integration is ever the first Brillouin zone . Defining the mean free path 


A(k) = X (k) v^cne has 

K = JL — ( d^k ^ cos‘^9 

(2 TT )2 


(2.47) 


and A- (k) is given by 
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i ^ 


2 TT Nq \ (1 - cos 0) (TpC©) sin 0 d 0 (2.48) 


where N,, is the concentration (small) of the defects, is the 

differential cross secticnat polar angle 0. The weightage (l-cos ©) 
measures the relative change in the component of the momentum along the 


(2.49) 


initial direction of motion. The quantity 

g- = 2tT (1 -COS 0 ) CJ^(©) sin 0 d 9 

is sometimes called the momentum transfer cross section. If is 

ipa,ependent of 0 as is the case for scattering from isotopic imputities, 
when the scattering is entirely S wave type, then ’CT'n equals the total 
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cr'Rs section which can also be cTOluated ty using optical theorem. Of 
special importance is the Hong wave behavior ofy\,(k). Since tJ* 2 (©) 
contains resonance denominators ^the quantity _/\_(k) registers a .gmall 
vaHue at the resonance in turn implies a smaller value of 

the the.ral conductivity K than when resonance is absent. Thus the 
reS'"'nance in (j p(©) implies a ’dip' in the thermal conductivity curve, 

In interesting feature of equation (2.48) is the interference 
be^'ieen the terms of different syiumetries . In particular terms of different 
inversion symmetries, (or parities) namely the ’gerade’ and the ’ungerade’ 
^erms show interference. This is seen veiy easily in the limit of long 
waves where one has for example f(k,l^) = f 4 - f cos 0 with f and f 
independent of 9. Then CTp = f* (k,!^^) f (k,^) 



whence by equation (2.49) 




4Tr 


g u 


■M- 

f f ) 

U g' 


(2.50) 


he oecond term on the right hand side shows the above mentioned interference* 


Thas is shown more explicitly in the next chapter- It demonstrates the 
imp-rfance of the weightage factor (l-cos 9) in equation (2.48). This 
point has also been discussed by Krumhansl and Matthew^"^. 


Equation (2.48) shows another important result that for .gtman 
concentrations for which the multiple scatterings from the defects can be 
neg-ected, the result for defects is just times the result for a 
single defect. If one starts from the consideration of random distribution 
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of defects at small concefnt rations and averages out the various 
quantities, one arrives at this ’'result as has been discussed by Klein 
and Vashishta^^. Thus the single impurity problem which appears too 
academic to begin with, is relevant to the pl^rsical situation under 
certain circumstances. 



CHAPTER III 


SEMPLE CUBIC UTTICE 


^ 3.1, The Model i 

In this chapter we consider a specific model and discuss the 

specific results in an explicit form. This model of a monoatomic simple 

cubic lattice with equal central and noncentral harmonic force constants 

4 

has been discussed by Montroll and Potts . There aiB N atoms of mass M 
each. The harmonic force constant characterizes both the central and 
non-central forces between the nearest neighbours only. This equality of 
central and non-central .forces enables one to discuss the problem without 
bringing in the question of polarization. With lattice constant a and 
Cartesian unit vectors e^, the lattice sites have co-ordinates 


R 


t 


a e ^ + m a e 2 ^ ^ 2-3 


(3.1) 


where , m and n are integers. We shall use either R or ( ^ ,m,n) to 

specify a lattice site. The cartesian components of displacement from 

equilibrium of the atom at site (L >m,n) are given by x 

4 

2 ^mn* discussed by Montroll and Potts , one can think of any of 

these Xi for example as a parameter associated id.th each lattice point 
Cmn 

( (,,m,n) and then this model is equivalent to a scaler model. With 
this it is possible to write the ecjiations of motion of the perfect 
lattice in the form of equation (2,1) 

- Mofe = 0 (3.1) 


/\ 

where A has the matrix elenBnts giTen by 


» *28 * # 
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imn, |^6 Q S in,m' S n,n' 

S m,m' S n^n'^S'l+l, C ^i-^> (£,11' |^^m+l,m' 

^ m-l, ^ m,m' ^^n+l,n’ ^n-1, n’ ^ (3.2) 


and 


JJ (k, imn) “ 


mn 


The same matrix element A(^mn^ Cm’ n’) will describe the motion when 

U (k, ( mn) are chosen to be -T/ or 3/ 

“ Lm *-i 


-mn 


The plane wave solutions corre^onding to equation (2.2) are 
(k^ imn) = exp |^i (k^L + k^m + k^n) J (3.3) 


which correspond to the dispersion law 


£0' (k) = -1^ (3 - s kp 


(3,4) 


Now let ttie defect atom be situated at the origin of the ^•oordinate 
system. It has a mass M and is coupled to its nearest neighbours ’ 

through force constant Y + AV* The perturbation matrix P has 
non-vanishing elements only when the indices refer to the impurity site and 
its nearest neighbours. D’=‘noting the nearest neighbour vectors by 

A 

the non-vanishing elements of P are 

P (2.S) = AM J - 6 A V , P (^, R ) = - ^ V 


and 


p (2 ,Eb> " *’ (fe'2) = A y 


( 3 . 5 ) 
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Here can ha^re compocents (+lj 0,0), (O^+l, O) and (0,0,+l). 

3.2. The Partial wave analysisi 

Is described in Appendix III, the reducible representation of the 

cubic group in the space of the affected lattice sites ( P -space) 
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reduces in the usual notation to 


n" « 2A, + F + E 

'5 Ig lu g 


(3.6) 


A 


This leads to the unitaiy matrix V given explicitly in Appendix . 

s 

III. Now using the procedure outlined in the last chapter one can find 
the scattering anplitude. Here we follow the notation for the Green 
function integrals described in Appendix I for the simple cubic lattice. 

Also the matrices ^p and ^ are explicitly written with same indices for rows 

A 

and columns as ere used to label the rows and columns of in Appendix HI. 
With these the elements of block diagonal form of (l+P^ , namely the 


tj, matrices of equation (2.28) are given by 


A 

t 




A A 





i/a / 

AV ° ° 



1/7. 


(3.7) 


where 


A 


H 


Q 


^ Aii- cD - 6 

AY 


t= ae + 

o 


^4l(l,l,0,Eg) + I(2,0,0,S3)^ 


- C +-L.J' 4I(1,1,0,E ) + 1(2, 0,0, Eg)] (3^) 


f p31 # # 


i I (l,0,0,Eg) 
^ I(0,0,0,E ) 


and 


and 


q[i - ^V{(aC^ 6eJ| 

■- h[(ve„)] 

(3,9) 

1 + ^ Si (0,0, o,s ) . 

- I (8,0,0,E )7 

(3.10) 

^Y 1 

J 


1 + ^ (0,0,0,E ) • 

+ I (2,0,0,E„) -21(1,1,0, Ejl 

C 


J 


Using these one can evaluate the scattering amplitude as discussed in 
the last chapter and the result is that the scatteilng airplitude 
f(k,ko) can be written as the sum of three teims 

f(k,k^) - fs(k,k^) + fp(k.4) + ^^(k^y (3-12) 

where the S-like wave amplitude is 

T I 

"s 




go(ko) 


+ l(2^0j0,B ) + l(0j0^0jE_) 


-6 AY + ^4 1(1, 1,0, Eg) 

)J. + 2 AY 1^1- l{l,0,0,Eg)J 

y ^ (=0= k. + C06 I(0,0,0,E^) j 

^oi 3 ] (3.13) 


The P~like wave amplitude is 


fp (k,kg) «= 


-?s^(.Ay) sin k^ sin. k^^^ (3 J.4) 



m «i33 m « 


and tbs D-like wave amplitude is 


DO 3 




D 


3(cos kg, - cos lc^)(cos k^g-cos ^ 03 ) 


+ (2cos k^~ cos kg - cos ^ 2)(2 cos cos k^g - cos k^^) 


(3.15) 


When the substitutional impurity is an isotope, the ? and D- wave 
amplitudes vanish and the S-wave amplitude remains in the form 




Sgiko) Amo 


Si 


where F„. = 1 - _ I(0,0,0,E_) 

2 Y 

This expression has been obtained by Callaway 


(3.16) 


(3.17) 


8 


3.3. The scattering cross section, its long wavelength limit 
and re sonance s : 


Here k specifies the incident direction. Ab already discussed, 
along the ^pmetiy axes of the crystal, coincides with the direction of 
scattering due to the behaviour of the constant frequency surface along 
these directions. This facilitates the use of optical theorem to evaluate 
the total cross section when the phonon is incident along one of the 
^mmetiy axes. Jllso in the limit of long waves, this can be done for 
ary incident direction. The scattering angle in this case is that 

between k and k . 

~ — o 

To demonstrate the use of optical theorem, we evaluate the total 
scattering cross section {5^ for phonons incident along (1,0, 0) axis 
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using equation (2.38). Then k = (k,0,0) and we set ^ = (k^jOjO) to 

evaluate the forwart scattering amplitude, k and k are equal in 

“o 

magnitude. With k^ = (k,0,0), gg(ko) = scattering cross 

2 

section can be expressed in units of a , where a is the lattice constant. 


V sin kL I Fg 


( Im D ReT^s - Re D Im Pg) 

V 1_AV 


(Im Fp) sin^ k - - J^A 0 (cos k-l) (Im F ) (3.18) 

^ I I ’^1 


where 


D = AM 60^ -6 AV+ 4 M of ^4 I(1,1,0,E„) + I(2,0,0,E^) 

2 V 1 


+ I(0,0,0,E^)J + 4 l(l,0,0,Eg) j. 

_ 4 i 1 - 1(0,0, 0,Eg) ? (2+ cos k) (2+ cos k) (3,19) 

A 2V 


^ ^ ^ (3-Eg) S(0,0,0,Eg) + ^73(0,0,0, Eg) 

+ 12 C(1,0,0,E ) - 4 S(1,1,0,E ) - S(2,0,0,E^)2+ Am»(3-E ) 

® 2V M 

X je ^C^(l,0,0,Bg) - S^(l,0,0,Eg^ G(0,0,0,Eg)^4C(l,l,0,Eg) 

+ C(2,0,0,Eg) - C(0,0,0,Eg)J + S(0,0,0,E^) £4S(l,l,0,Eg) 


+ S(2,0,0,E^) - S(0,0,0,E^) 


J] 


(3.20) 
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Im (S-Eg) C(0,0,0,Eg) + 


jQcf... ? 


M 


2y L 


70(0,0,0, Eg) - 123(1,0,0, Eg) 


- 4C(l,l,0,Eg) - C (2,0,0,E„)J + 


A 


- 3 / p. ■ ^ ^ (3-Eg) 


2V 


M 


Ee = 1 + 

2 2 


X ^-120(1, 0,0,Eg) S(1,0,0,E3) + C(0,0,0,Es) ^^S(l,l,0,E3) 
+ 3(2, 0,0, Eg) _ S(0,0,0,Eg) j + S(0,0,0,Eg 
+ C(2,0,0,Eg) - C(0,0,0,Eg) j J 

S(0,0,0,Eg) + 3(2,0,0,Eg) J 


(3.21) 

(3.22) 


In 


|^^C(0,0,0,Eg) + C(2,0,0,Eg) 


(3.23) 


Be Fp = 1 + -^-^^3(0, 0,0, Eg) + 23(1, 1,0, Eg) _ 3(2, 0,0,Eg)J (3.24) 

3 


Im F, 


= -^|^C(0,0,0,Eg) + 20(1,1, 0,Eg) - C(2,0,0,Eg) ( (3.25) 


Here we have used the C(p,q,r,Eg) and S(p,q, r,E ) defined in Appendix I 

28,32 ’ 

and their values have been tabulated. 

In the special case when the substitutional inpurliy is an 
isotope, AV= 0 so that D = and = 1 - ^ l(0,0,0,Eg) 

which leads to the result 

We define 0) 


■ 1 f 6 k \ o)^ ^ . 

V sin 4. 2V y jp 12 (0,0,0, Eg)) (3.26) 

2 2 V I Sij 

(3.27) 


oS 


M 


and in the limit of long waves 


^ I. Jl(0,0,0,E,)j = ^ CC0,O,O,K,) . 


and then with sin kH^k, equation (3.26) yields 



,. 35 ., 


O', 


1 /Am s2 0£>^ 

— — tk / ~~ y 

l\ 


M 




in units of a. (3,29) 


This is exactly the expression obtained by Callaway . The scattering is 
purely S-wave type . 

The long wavelength limit of (Tj equation (3.18) is given by 




+ 






k‘^(lm F^) + 


Im D Be Fg - He D Im Fs 


] 


aV 




(-Im ^) 


(3.30) 


where the quantity D now becomes 


D = -(k^ /6) +^i4MJ^C4l(l,l,0,Eg) + l(2,0,0,Sg) 

2 V L 


i 

v2 2 

+ 71(0, 0,0, Eg) - 121(1,0, 0,Eg)J+ ^ ^1(1, 0,0, Eg) 

- 1(0,0, 0,Bg) ^ (3.31) 

where l(p^qjr,Eg) are to be evaluated in the limit of long waves* 

4 

Retaining terms upto the order of W ^ the D-wave amplitude does not 

6 

contribute because it is of the order of C<^ and hence 

.2 . ^ ^^-,2 


w |Fgf 0 '^ i J 


A)oS^ 


(3.32) 


in 


in 


units of a^ where the constant b This result shows that 

P 

the long wavelength (low frequency) limitj^ the scattering cross section 
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is proportional to oj showing Rayleigh type of scattering, provided 

that there is no resonance which is seen hy examining Re Fg and Re Fp . 

The S-wave part of the cross section is affected ty both /i M and AV 

whereas the P-wave part depends only upon A"/* 

It is worth while to point out that the number of partial waves 

contributing to the total scattering cross-section defends strongly upon 

the direction of incidence. This is discussed in some detail in the 

next chapter. For example if we evaluate the scattering cro.-vs section 

for the ( 1 , 1 , 1 ) direction of incidence the D-wave part will not contribute 

to an expression like that given in equation (3.18), 

A resonance in the scattering cross section will occur for that 

value of the frequency at which any of the quantities ReF 3 , Re Fp or 

<= is small and positive, 

(Ka F)' 

The real and imaginary parts of Fg, Fp and F^j are given by equations 
( 3 . 20 ) to ( 3 , 25 ) .f One can study the resonances by taking given values 
of ^ M and ^ and examining the criteria for resonance. The 

M -y 

figures 1 and 2 give the results of such a study. In these the value 
of the parameter E , defined in Appendix I is found at resonance for given 


values of 


A M 


and ^ Y . We find that the ? and D resonances depend 
upon 4 » only. For sufficient decrease in force constant the P and D 
resonances appear at the low frequency side of the band. If this is 
increased the resonances appear towards the high frequency side of the 
band. Indeed in figure 2 if one changes the sign of ^3^ 
corresponding Egj^ simply reverses the sign having its magnitude intact. 

The S-wave resonances depend upon Am as well as AV . For a 
large increase in mass, resonances appear towards the low frequency 
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(large positive E^) side of the hand, no matter what the change in 

force constant is. On the other hand, for small changes in mass, the 

change in force constant plays a dominant role in determining. the resonances. 

4 

At resonances, the scattering cross section does not have a 

dependence in the limit of long waves. Since the real parts of Fg and 

Fp become zero and the imaginaTy parts of Fg and Fp are. proportional to 

f y * then equation ( 3 .30) for 0""* shows that it is proportional 

^ . .-2 

f ^ a't resonance. Further the width of the resonance is 
proportional to J • Thus the value of TOaomnee becomes 

laige .iy a factor with small Ci) • Using the expressions given 

by equations (3,20), (3.22) and (3.24), ojfe can study the localized modes 
of S,P and F symmetries. 

3,4 Phonon scattering and the mean free path ; 

We evaluate the mean free path of phonons in the limit of long 
waves by usihg equation (2,48) . The differential sjiettering -eaoas eectioj* 
in this limit is simply tho s<jiia-re of the scattering ajj^ilitudo 

f 

which is given by 



with b fid 1.14 and Q is tho scattering angle. 

Then 

T 

' ,2 

I f(©) I (1-c.os 0 ) sin Ode 

0 
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The last term on the right hand side is tho interesting tnim showing 
tho intorf'^ronco between terms of opposite 'parities'. Hero 'paridy' 
signifies the behavior under inversion at the centre of inversion. The 
S-type of tem does not change sign under inversion and it is the so 
called 'gerade' teira and the ?-type is'ungerade’. 

At resonance the cd ^ proportionality of -A- does not 

_2 

hold and the dominant term can be shown to be proportional to Q) 
as in the case of (5”^* 

Thus the simple model considered here shows the important 
features of the scattering of phonons except the question of polarizaticai 
viiich is discussed in chapter V, The long wave length Rayleigh "type 
of scattering and the modifications at resonances are explicitly 
obtained here. The directional dependence of the scattering cross 
section, the interference of terms of different symmetries in the meanfiee 
path expression are the other important features of the ]»rocess and 
these are further elaborated in chapters 17 and V, In the following 
section we consider a model of diatomic simple cubic lattice to 
demonstrate the ease \d.th which the results for the mono-atomic lattices 
can be cariled over to the study of diatomic lattice. In fact this feature 
adds to the usefulness of the results for mono-atomic lattices. 
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3,5 Diatomic Simple cubic lattice 

Let a diatomic simple cubic lattice be consideied viiose 
atoms of type I have mass each and those of type II have mass 
]y^ each. The central and non-central harmonic force constant 
between nearest neighbours is denoted by Y . One can write down 
the equations of motion and perform the so called M transformation 
used by Montroll and Potts^~® so that the equations of motion for 
the diatomic lattice becomes identical to those of monoatomic 
lattice with M q) in place of M vhere, 

M* = 6 T + >/ (M^ C()^- 6Y) - 6i) (3.35 ) 

Then the dispersion becomes 

M* £0^(k) = ey - 2y 2 cos (3.36 ) 


. cos k. 
i X 


and this leads to 

s/(M^cO^(k) - 6 y ) (M2£«)^(k ) - 6 y ) - - 2y 2^ 

(3.37 ) 

Solving this equation for ^(k ) one gets the optical 
and acoustic branches as, _____ 




+ -Jx. ) f 1 + 

^ ^ L J 


and 

O)^ (k) -= 3 1 ( 
ac 


-L. + — i_ ) 1 „ 


4 ( 1 - i cos k^ ) “^ 

* -k- - 

. (3.38a) 


4 (l « -i- cos k. ) 

^ ^ — ^ ^ 

- J 


(3.38b) 


The maximum frequency in the optical band is Ct>^ 


ipM 


i 1 

® y minlinum frequency in the acoustic band is 


# # 4:2 • # 


2 

cO = 0 . The miniinain in the optical band cD and the 
a cm opm 

maximum in the accufitic band decend upon the relative 


magnitude of and Mg . ^or > Mg , = -2j_ and 


. 6 

whereas for M^^ M- = — — and c5 ^ ~ -n • The gap 

!» opm M^ acM Mp 

between the two bands is of magnitude 6 'j^|( 1/M^- 1/M^)|. 


In analogy with the monoatomic case the Green function integral 


becomes , 


■j , S ' exp 5 i(lc, t + k m + k n ) Z 

G ( cPi L m , n ) = - 1/N ^ II I 1— ' 

i H cO (K ) - - 1 a 

^ i (k^ t + kgffi + k^n ) J 


2^ 


exp 


^ V ^ k E - 21 cos k. 

i 1 

9r 

where E = - ~ J (M^ Od^ - 6 ){U^cO^ - 6 yT 
1 


= - 3 



(3.39 ) 


(3.40 ) 


(3.41 ) 


C02 _ 
^ ' ojs" 


for M^ = Mg . In terms of the E defined 


os 


It may be noted that at each step one can check the results with 
that for the monoatomic lattice by setting ^ ~ ^2 * quantity E 

goes over to E 

s 

above , 3 ^ E ^ 0 defines the acoustic band and 0^ 'E'^ - Z 

defines the optical band . Within the gap, E has an imaginaiy value. 

These have been discussed in detail by Mitani and Takeno 

OJqpm 
2 



u-J opm 
2 

cOacM 

2 


6^: 


acm 
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Thus taking aopropriate values of E one can consider the 
scattering of phonons in both the acoustic and optical bands. Optical 
phonons do not contribute to the low temperature thermal conductivity, 
but it is interesting to consider the scattering of optical ohonons from 
the substitutional impurity. 

Consider a su'hstitutional impurity put at a site of an atom of 
type I x'hich may be taken to be the origin of the coordinate system. 

Let its mass be + A M and the force constant be "Y 

The elements of the perturbation matrix after the transformation 

get modified into the following. 


P (0,0 )=^( MQ)^-6Y), p(Q,R^)=p (R^,0 ) = Ay. 


P (R , R ) = - 


hy 


Following the method described in this chapter ^ after oerforming the 
M transformation) one gets the following resona'»^ce denominators also 
obtained hy Mitani and Takeno^^. 

Fg = 1 -H - ( AX _ G ( a)^, 0,0,0 ) 

+( ) ( 1 + AL ) M ci)^G(c;>^l,0,0 ) 

V 1 


(3.42) 


"1 


and Ftj = 


- ^G ( 0)^,0, 0,0 ) - G ( a>^,2,0,0 (3,43 ) 

1 “ 0c)^»0,0,0) + G ( 0)12,0,0) J. 


- 2G(60^, 1,1,0) 


where 5 - 6 7 ) / (M^ Cd^' - 6 y ) 


(3.44) 
(3.45 ) 


bithin the bands, the parameter E is real and the Green functions 
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o^curing above can be written in terms of I ( P,q, r, E ) as, 


C ( p,q,r,) 


■27 


I ( o,q,r,E ) 


( 3.46 ) 


which in turn can be written in terms of the tabulated functions 
^ (P>Q>r,E) and S (q,q, r,E) . Here cD * is related to E by the relations, 
(i) Optical band, 

.2 


! oom acM ^ \l opm acM ^ 


+ 4 ”1 

opm acM 


and (ii) Acoustic Band, 
2 


cb 


i [ (cd^ + cd^ - I 

‘ 1_ ^ opm acM / 


. ^ ^ Cd'\''. - I { - O)^ )' 

opm acM ^ opm acM 


+ 4 iE^/9) cb^ cd^ 

opm acM 


] 


( 3.47 ) 


( 3.48 ) 


Using the above expressions for the resonance denominators 
hxtani and Takeno^^ have investigated the localized modes both above 
the top of the optical band and within the gao between the acoustic 
and optical bands. We have investigated the resonances in the acoustic 
and optical bands and the results are given in the figures 3 to 6 for 
the cases; (i) , p = M^/M^ = 2 and (ii) , p = M^/M^ 

= -i- . These are the same as those used by Mitani and Takeno^^ for 
investigating the localized modes. 

The change in mass does not affect the resonances in ? and 
D like waves. Resonances only in the optical band are shown for 
these waves. The result is veiy similar to that of the monoatomic 
lattice case and as the force constant increases, the resonance 
frequency moves to-wards the high frequency edge of the band. 
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Figure 3. P and D-wave resonances for Mj> and p- M 2 /Mj_i in Optical band for 
diatomic simple cubic lattice. 


i 
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Figure 6. Typical S-wave Resonances in the Optical band of a diatomic simple cubic lattice, 
lattice. Curves a, b, c and d are for M|> M 2 , p= 1/2 with Ay-y=-- ^ 

0 and 4*, 1 respectively and curves a‘, b’, c* and d' are for ^ 2 ^ ^1* with 

/^y/y=-,5, -.3, 0 and -t. 1 respectively. 
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In all the figures E denotes the value of E at resonance. For 

E 

the cfse M / M which implies that the substitutional impurity has 

X ^ 

been put at the site of an atom of mass (heaTier mass), comparatively 
low values of increase in force constant proc^uce P and D wave resona- 
nces in the optical band. For higher values of increase in the force 
constant one gets these resonances if an atom of mass is substituted 
by the impurity . 

The S - wave resonances are interesting. Tr the acoustic band 
one gets the usual result that as the increase in mass becomes large, 

one gets low frequency resonances, no matter what the value of the change 

in force constant is. This remains true for both the cases of (i) Mg 
and (ii) M M . But in the optical band one gets different behaviors 

for these two cases. the site of the heavier atom is substituted 

Mg) f one gets resonances for a decrease in mass and the smaller 
iiia.sses give rise to resonances nearer to the bottom of the optical band, 
^^en the site of the lighter atom is substituted, one gets resonances 
for which the smaller masses give rise to resonances near the top of the 
ootical band. It is interesting to note that from E=OtoE=-l 
large decrease in force constant is required to get resonances ( not 
shown in the figure). This can peihaps be attributed to the behavior 
of tho density of phonon states in this interval. 

From the above analysis it becomes clear that while in the 


monoatomic case, only heavy masses can produces resonances and light 
masses give rise to localized modes, in the diatomic lattice light masses 
can n reduce localized modes above the optical band and resonances in 
the optical band. In addition one can get localized modes in the gap. 

f 

43 

The present analysis, along xd-bh that of tfi-tani and Takeno leads to a 
rather complete analysis of the problem for diatomic sinple cubic lattice. 



CHAPTER rV 

BODY CENTRED CUHTC AND FACE CENTRED CUBIC LATT. GES 



j this chapter the results of a study of fhe proble:i of 

scattering of phonons from a substitutional inpuiity in a siEolified 
mo'^el of face centiBf^ cubic anc* body centrer^ cubic lattices are 
presented. The raothod of analysis is the same as that presented in 
chapter II. The displacements of the atoms in the lattice are assumed 
to be scalar parameters Tdiich depend upon lattice sites. This is not 
realistic for lattice d.yn£.inics as there would be no polarisabion of 
phonons in this model. However, some essential features of the process 
of scattering are brought out oven in this model end one oan get a 
qualitative insight into the oroblem of scattering in more realistic 
lattices. As shown in chapter T, the results of the present chapter 
can be obtained by a suitable choice of force constants in a realistic 
model. The present analysis is, however, oomPlotoly applicable in the 

problem of scattering of spin waves or band electrons f -om .substitutional 
i’ impurities. 

In the two lattice models considered here each atom has 
mass M and nearest neighbour harmonic force constant Y" , except 
for the impurity atom which has mass M + AM and is coupled to its 
nearest neighbours through a force constant Y- AY . The impurity 
is taken to be at the origin end the lattice vectors R have 
Cartesian conponents^^, 

f! Sb *= (t + m/2)a £2 + ( ir/2 )a + (n + ^^)a (4.1a) 

I, 


• # 
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r.r.'d 


Sf 'tCL'*' Bi)a ©o "*" TCni '•' n.)a * ( 4.1b) 


be IB c is the length of the side of the unit cell and e£^ are the 
unit vectors along the Cartesian axes, and and n are integers. 

The total number of atoms is N. The subscripts 3 and ? refer to 
Eody Centred Cubic and face centred cubic lattice', ic Elective ly. 

The fiequencies (squared) are given by^ 

2 2 ^ ^3 

® i*J (l - cos -4— cos — ^ cos ) , (4.2a) 

B OB 2 2 2 ^ 


and 


2 

P 


k 

2 


Co (k)'=*4; (3~ cos — cos • — cos — - cos 

■"OF 


k. 

£ 

2 


k 

£ 

2 


2 


h \ 

— cos cos — =— ) , ( 4,2b) 

2 2 ' 

vhere to *=8 "//M, tO =4 ~f /M, and k. are the Cartesian 
OB OP 1 

components of the wave number k. These can be obtained by writing 

the equations of motion in the form of equation (2.1 ) with the 

A 

components of A given by , 

B,F 

^ fm* ,nd-l ^mt,nwL' 

" S 1[+1 2 m» ;m-l ^ , 1-1 ^ m^m+l ^ 


C ft ^n^n+1 ? n> ,n!~i ^ 

/ ^ 

^L[L+i f m*,in«l ^n*yn+l 

» f -1 ^ m’^m+l ^ n’ ,nr4 ^ ^ 

S7ttl 


(4.3a ) 
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Imnj, t'm'nO = yri2 C (T - f ^ jT , (£ , +C ) 

^ ‘^{‘( dm’m 5n»n (o m'm^ ^n’n+1 ^n'n~l^ 

) y l> ? ; 

f Sn’n "£ m'mfn'n^^M+l 

J ^ V J J ^ j 

^ ^mj^iir+-l*b i^n-1^ ^i^Et-l ^njn+l^'“ (+lC rijn-l 

‘^^(.'(-1 S'n^n+1 ^ 


^ n^n ^ ^pj+l£ m^m-l ^ p [~>1 S'mf m+l ^ J 


(4.3 b) 


The pertuibation matrix P of equation (2,5 ) for these cases 
he.s non— Vanish.! 3ng dements only when the inciices refer to the impurity 
site and its nearest neighbours. ¥o shall donote the nearest neighbour 
vectors by ^ The non— vanishing elements for the two lattices are, 


p (£,s) = 4/ , p (R^,^) .-^Y> 

P (2.S„) -P (S„,2) = , c 4.4 ) 

where Z_ = 8, and Z =12. 

O p 

The unperturbed nonaal modes are given by, 

"= IT^ exp (i k.2/e )' ( 4.5 ) 

The elements of the Green function matrix, along with 
their low frequency expansions are given in Appendix I, 

I .4: ,2 The Partial wave analysis 

Here the point group of the lattice is the cubic group (^. 

The representation of Oj^ genereted by the displacements of the 
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affected lattice sites reduces, in the usual notation to, 


n „ = 2 A, „ + + A , 

' B Ig lu 2g 2u 


(4.6 a) 


and 


= 2 A, + F^ + E + F^ + F^ 

* F Ig lu g 2g 2u 


(4.6 b) 


One can use the method of projection opeiators to obtain the 

A 

symmetry vectors which constitute the columns of V matrix as 

A A 

desciibed in Appendix III, where V and V are found explicitly. 

B F ^ 

The indices of rows end columns of the V matrices as 

/\ A 

specified in Appendix III for V and V are used similarly to 

. * B F 

A A A 

write the matrices Pg, Pj, and in their explicit forms. 

The inverse block matrices t^ occuring in equation (2.28) 
are as follows. For the body centred cubic lattice. 


\b “ F 


SB 


A 


1 + A M to 

- giB> 1 ^ ^oB ® ^ hB 


AAA 
^2B " "*^38 "" “ "F, 


A 


A 


and 


PB 


, t *= t *= t = 
" 5B 6B 7B 


A 

t 


8B 


“^FB 


BB 


(4.7 ) 


Here, 


Ag *= ■“ B a/ . 


”b ■= Ab SlB * ^7" (3 «2J * Sgjg + gQB hs") ’ 


(4.8 ) 
(4.9 ) 


^SB - - ”b) (1 * Sob ^ 3 AV ) 

- 8 mj a/ (g^B - %B> ' 


(4ao) 
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^PB “ ^ ^ ^2B " Sb ' ^4B ^ ’ 


and^ 


- 1 


./( 


^DB ~ ^ “ ^ “ ^on ~ ^-TTS "*■ S^-D ) > 


OB 2B '='3B 4B 


- &■/(( 


■FB - - “ » v^OB " ^^2B ^®3B "^46^ * 


(4.11) 

(4.12) 

(4.13) 


The notr?.tion for the Green functions g is (iescribef^ in 

iB 

AppencJix I. With these^ the scattering amplitucile is, 

’ -o> ^SB • K') * fpB (i’ * Sb 


+ f^(k , k ) , 

jigv-. ^ > 


(4.14) 


who re , 


^SB 


^• 0 ( 1 ^) 


=B'^' j , 2 


•SB 


[ (1 + m tO“gj^g - Jlj) (Ig 


8 A/c^GjC,) C^C^Cj) ■. 8 0^,0„2C^3- Af X 

{(1 - * ^B^CB " ^ 'iB > 




j] 

[ ®d2"3®ol‘'oB'=o3 


(4^5 ) 


+ S C C.S „G _G , + S_G.C_S -C ,C ^ 

2 3 1 o2 o3 ol 3 1 2 o3 ol o2 


f (k, k ) = 


DB F, 


EB 


[=l®2S3'=al^o2So3 * 


+ c_s, s^c _s , s ^ , . 

3 1 2 o3 ol o2 f * 


] 


(4.17 ) 


and. 




(4.18 ) 
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Here 


C ■= cos (k./2), S. = sin (k./2), C . = cos (k ./2) anf’ 
1 1 01 01 


- sin (k /2). The expression for g_(k ) is given in Appenf^ix II. 
Similarly, for the face centred cubic lattice one obtains, 

2 


^SF 


1 +' 

- 12^Ay(gQp - g^p 1 + + 12Ayg. 


A 

A 

A 

^2F 

" "'^SF 

'^4F ’ 

A 

A 

A 

■^TF 

"^SF 

“ Sf 


K 


TF 


DF2 


"5F 


, and t 


A 

t 


6F 


F^OF 

1_ 
^DFl 



A A A 2 

lOF %1F ^12F'"T^“ 


where. 


Aj, = Am _ 12 A y , 

“F - VlF *AV(V * «SF " * V ’ 


(4.20) 

(4.21) 


= (1 + ZlM - H.J, ) (1 + Ayggj, + 12 AV^ip ) 


12 /iynip (Sqj, — ) 


F _ = 1 + A 7^ - g . - 2g + 2g + g ) , 


PF 

'DFl 


0^ IF 4F 5F 

^ ^ +Ay (- gQp + 2g^p - 2g2p + 2g^p - gj-p ) , 


and , 


F =1 +Ay(- g + 2g - g ) 
DF2 ^ ^ ^OF ^2F ^5F ^ 


^Ffi’ ^ ■*'Ay(“ gf.-Bt 2g^^ - 2g.* + g ) , 


FF V sqP ‘^&i5' ^4F “5F 

The scattering amplitude in this case can be written in the form, 

fj,(k ,k^ ) = fgj (k.k_^ ) ^ (k,k^ ) + (k , k^ ) 

" "dF2 " "ft '^o > ' 


(4.22 ) 
(4.23 ) 

(4.24) 

(4.25) 

(4.26 ) 


(4.27) 
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where, 


■f fir V ^ 

o-m 


(1 +AM u> _ mp) |Ap + 4A-f( C^Cg 

+ CgCj + ^3^1)^ 2ip(3 - ^1^2 ^2^3 *^3^1^ 

" ( =01 <=02 * ‘=o2<=o3 * h.3=ol ^ X 

[ 4 " * J'-F^OF * l^Aj’gjj,) {3 - - C^C^ - C^C^ ) 

(Sgp •“ ^ ^2^3 ^3^1 (4*28 ) 

Sh‘=2 ^ '^3)Sol(Co2 + <=03) 


fpF(i -!£o) = - ^ %(k„) 


+ 82(03 + Cp)Sq2(‘^o 3 ^ol ) 

F S3 (C, + Cj CJ ] (4.29 ) 

f^n^iE-fc) " -f ^ r=(Sh - h'=2X'=o3'=ol- =ol'^c2> 

+(20203 - 030^ - OpOg) (20^20^3- 0 ^ 3 C^p- C^^C^)J ( 4 . 30 ) 

^DF2^^ .Mo) == gp(Mo) [SiVol®o2 ^ ^ 2 %^o 2 ^o 3 

and, 

>V--I~ h %^bl <=2 - " 3 )Soh=o 2 - <=03) 


+ 22(23 - Cj ISgj^E^oS ■ ‘^ol) 

" Ss'h - <^ 2 ^ 3 ^%! - ^ 2 ) ] • 

Here we have distinguished between the two D-wave 

amplitudes fjjpp and £^^2 arise due to the two irreducible 

representations E and of 0. . This identification is described 

g 2g h 
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in the character table give’^ in Ae-oen<^iy III. It nay be mtc'^ 
that above is veiy siinilar t^ for the simple cubic 

lattice given by equation (3.15). This may be expectc'^ '^ue to the 
fact that both of them arise fue to the irreducible representation 

^h ‘ is v-eiy similar to fg^^, both 

arising from the irreducible representation F , 

For the isotopic impurity case only the S-wsvc amplitude 
remains in both cases, and the form of the amplitude is the same 

as that in the simple cubic lattice given by equations (3,16) and 
(3.17) . 

^ scattering c r oss-section, its long ypvelength limit and 

resonances 

As iKsntioned in the last chanter, k specifies the incident 
direction and 1^, the d.irection of scattering whenever k is 
along any of the symmetry axes. The total sca+tering cross-section 
is found by using the Optical Theorem. An interesting feature of 
the scattering process is that the number of Partial waves that 
contribute to the total scattering cross-section depends on the 
direction of incidence of the ohonon. This is expected, in view 
of the fact that the ’potential’ due to the impurity is a non- 
central one, and its 'range' depends on the direction in which 
it is viewed. Ve will discuss the results for the total scattering 
cross-sections corresponding to a few incid.ent directions to 
demonstrate this effect. 

The forward scattering amolitude, which is necessary for 
the use of the Optical Theorem to evaluate the total scattering 
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cross S0cti-r, is epsily foun^ xjhen the f^irectic'^ of 3.r>ci'-''ence is 
(loo), (no) nr (ill). Cre can then set k = . For incirisnce 

along (lOO) direction, 

k = = (k, 0 , 0 ,) 

ani in the Body Centred Cubic lattice, 

1 


and 


gs^feo) - ~s^-y- (k/ 2 ). 


^j^(lOO) _ 1 

TB 4^ Sin k 




- (Ke Fgg) ( . 




(4.33) 

in units of p^. Here Dg^ is given by, 

2 

®B1 ^ iiMiO g^g -“b)(^B + 86 /cos k/2) -86'}Jng(l - cos k/2) 

+ 86V cos(k/2) ^(1 - cos k/2 )(1 + Agg^g + 86/g^g) 

- (Sq^ - g^g) (ig + 84ycos k/2 ) j* . (4.34 ) 

cr"(ioo) 

jg shows no D and F wave contributions. For (ill) direction, 

i 1 


k = kn = k/ 3 ‘, k/32 ), 


and , 


(ill) 


( 5 ^ ^ 1 

2y(l2)Vsin(k/l27) Cos2(k/l2^) 

^ ( Im Dg 2 )(lle Fgg) - (Re Dg^ )(lm ^sg)J- 
cos^(k/l2’^) sin^(k/l2'2') (in Fgg) 

- cos2(k/l2-^) sin^(k/l2-^) (Im F_ ) 

1 DBl 

sin^(k/L2'^) (imF^)"] , 

I ^FB p ^ J 


4 - 


(4.35 ) 
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in units of a . Here ^g 2 ^iven by, 


Db 2 = (1 + AM - mg^CAg + HaY* cos^ k/l2^) 

- 8 111^(1 — cos^ k/l2^ ) + Sii\ cos^ k/l2^ X 

[d - oo3=k/ls^) (1 + + e^iygjg ) 


k/l2" ) X 


(4,36 ) 


Thus rjt-Q ^ shows contributions from alb the partial waves. 

For (no) direction, it can be shown that the F wave 
contribution is suppressed. 

In the face centred cubic lattice the corresponding 
expressions for the total scattering cross section can be found 
in the same way. For example. 


<g^(iio) 1 

TF 


2Vt8*)^(sin k78i')(l To os k/S'J) 


”f "IF ^ 
SFl 


.^(Im Dp) ( Ee Fgp ) - ( Re Dp ) (im F^^ )^ 

sin^k/8® (l + cos k/8~)^(lm F^ ) 
t ^ irF 


cos^k/82 (1 _ cos k/P^)^ (im F ) 


\^?f\‘ 

- _L^_V 

- ---"C sinVs^ (Im r ) 

F DF2 ^ 

' DF2t 


DFl' 


^ sin^k/8^ (1 - cos k/82)^(lin Fp^ ) 


I F 

‘ FF* 

, . 2 

in units of a , with D^ given by, 


(4.37) 


JL 

Dp = (1 + (i^2g^p - mp) j’Ap + 4AV(2 cos k/S^ + coaVs"^)} 

— 4A'y®p(3 “• 2 cos k/8^ — cos^k/S^) + 4AV(2 cos k/S^+cos^k/S^) X 
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Cl i O i 

j — (l + ApgQp + 12^g^p(3 - 2 cos k/8- - cos^k/S') 

- (gQp - H 4^V(2 cos k/fi-^ + cos"^k/82)l S (4.38) 

Similar oxprossiors can be obraired for incidence along (lOO) and 
(ill) directxons . Table I gix^es the number of partial wares contri- 
butirifj to the total scattering cross section for various rlirecticns 
of incidence in simple cubic, body centred cubic and %ce centred 
cubic lattices. 


Table T 


Crystallograohic 
«xes of incidence 

Partial U^vos Contributing to C>t 

Simole cubic 

B.C.C. 

F.C.C. 

(100) 

S, ? and D 

S and P 

S. P and 

(110) 

S, ? and D 

S, ? and D 

S. 

(111) 

S and P 

S,?,D,F 

S, i end Dj 


In the limit of long waves the dependence of the total 
scattering cross section on the incident direction ceases to exist, as 
expected, and the exoressions for the two lattices are ( in units of 
a^), 

4 

ttK 




SB 


Oj 


op 


+ — 


16 


V 


<n 


TP : 5“ 

47r ^ 

I SF 


"PBr 


(4.39) 




SFI 

1 




QloF 


/ cD \4 / ^ y n2 
2“ ^ 


V 


•^OF 


t 


(4.40) 
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TItj. scau'-orin^ cross s-'ctrcn in both cases reveals Hayleigh— type 

scattering when there is no resonance. F— and aie to be 

bF L F 

evaluate'^ in the low ■freqiiGncy limit* The cliai^ge mass affects 
only the S wave part, and the change in force constants affects all 
the Partial waves. 


As discussed in the last chapter, this Ha 7 /leigii-typc of 

scattering does not hold at resonances. At a certain frequency 

the total cross section will be prooortional to CO in absence 

r 

of resonance but at resonance it will be much larger and xd.ll bo 

p 

proportional to C(J ~ . This is because xjhen there is a S-type of 

r 

^resonance, the Ee g p will vanish and Im Psg ^ is proportional 


to 


f _il_ 3 
h cb OB,F J 


so that ^ ^ is 

Kj TB,F 


pr'~''Dortional tof ----- 1 

^ OB,/ 

Thus the total cross section becomes laiger by a factor / 

^ CO 0B,F/ 


-2 


The same is true for the ?— wave resonances where Re F^ 


0 


and Im F^- is prooortional to / 

*<- CO OB,F>' 


The width of 


these resonances will be gi-en by equation (2.43) . It is seen 

easily that fj-, is prooortional tn / — ^ , 

V CO 0B,F/ 

The resonances in the scattering cross section can be 
investigated in the same way as in chapter Til. The figures 7 & 8 
show the results for body centred cubic lattice. Here ELg denotes 
the value of Eg ( defined in Appendix I ) at resonance . 


From the figures one can see that the , D an:! F- like 


resonances occur only when the force constants are decreased 
sufficiently. These of course do not involve the change in 
mass. The S-wave resonances depend on AM and iiV- For 
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slight increase in mass and sufficient decrsa.s''- in the force 

constant ''ne obtains resonances. For high value of the 

increase in mass, res'^nancos eopear near the band edge, no matter 

whether the force constant is increased of decrease'^ whereas for 

slight changes in mass thr change in fcnace constant plays a 

dominant role. The localized mod^s can be investigate'^ using these 

resonance (^encminators outside the ban'^ • 

^ ^*«'4Thonon Scattering and the mean free Path 

The mean free path of phonons can be compute'^ using 

equation (2.48) . This involves (5^ ( 0 ) and the scattering 

angle G . OT'p (O ) is readily computed from the scattering 

amplitudes. In the long wavelength lirit only the S and, ? terms are 

4 x-t 

significant upto the order Ct> • The expressions for the 
scattering amolitudes for the two lattices in this limit ere as 
follows . 


(a) Body Centred Cubic lattice: 

' <"> -[4- B© {-fr) (‘ 


• ■“) - 


y 


ir VFp3 ; 


At 

V 


cos 0 




(4.41) 


which, when substituted in equation (2.48) , gives us. 


—i— = 16 N a^ 


37r]Fpij: 


(vo^b) ~7y 

(-vf -A 
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i SB 

(b) Face Centred Cubic lattice: 
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ant^ Drocee'^ing in the sane way we obtain, 

-L = F a 2 (^J±Y r — I 5- fAJLf (i-, J: 

A(cb) ° !-47r|FgpjS L My k Y j 37r[F^j,|2 
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(4.44) 


It ma"'" be note6 that in the expression for the mean free path 
the two partial wayes lead tn an interference term like the one 
grton by equation (3.34) for the simple cubic lattice. 

In conclusion, it nay be stated that even in the somewhat 
unrealistic models of cubic lattices studied here, some important 
features of the scattering process of phonons from a substitutional 
impurity becoire evident. The total scattering cross section depends 
on the direction of incidence of the phonons, as also the number of 
Partial waves that contri-bute to the cross section. But in the limit 
of long waves, no such dependence is there due to the isotropy of 
the constant frequency surfaces. The dispei*sion laws of the models 
considered here lead to Debye model in this limit. The mean free 
path for long wave phonons contains terms that arise due to the 
interference of S and ? partial waves. A comparison of the results 
obtained in this chapter with those obtained in the last chapter 
shows a great deal of simile rliy among the three types of cubic 
lattices. 


CHAPTER Y 


A MOREL 'aJITH PHONON POLABIZATION 

r~ 

'? 5 .1 Introduction 

The proceee of the scat+ering of ohonons from a substitutional 

impurity considerpc! in this chapter includes the polarization of 

phonons in the lattice . This problem has been discussed briefly by 

45 

Thoma and Ludwig for the isotope defect. Maradudir has treated 
the isotope defect case in some detail and derived the scattering cross 
section in the limit of long waves by using Debye approximation. Here 
we treat this isotooe defect case in detail without using Debye 
approximation in the limit of long waves and then we discuss the 
problem for a general substitutional impurity with grouo theoretical 
partial wave analysis which has been done for the scalar models -in the 
last two chapters. 

The problem is mathematically involved in the senee that 
significant part of the calculation requires cumbersome numerical work. 
Physically the process is complex because here the concept of polarization 
of phonons is not as simple as in the familiar case of photons. In 
general each state of polarization of phonons corresponds to a 
dispersion law peculiar to that state, unlike the case of photons 
where all the polarization states correspond to the same dispersion. 

All these dispersions become identical if for example the model is 
such that one Cartesian component of the displacement of aW atom 
is not influenced by the other Cartesian components of displacements 
of the same atom or its nei^bours. In chapter III we considered 
such a model. 
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The scattering problem can be stated as follows. Let phonons 
with a definite polarization and hence corresponding to a definite 
frequency branch be incident on a defect site and get scattered. In 
general they will get scattered into other states of polarization 
also. As expected the scattering process depends strongly upon the 
incident direction. It is requir^^d to find out the scattering cross 
sections for the process. In the model that will be consider'^d here, 
the constant frequency surface at low frequencies will not be spherical 
as is the case with the usual assainption of Debye aporoximation in 
this limit . 

S 5 .2 The model 

A monoatomic body centred cubic lattice is considered in which 

each atom has mass M and is coupled to its nearest neighbours 

through central and non-central harmonic force constants ^ ^ and g 

46 

respectively. Such a model has been discussed by Fine and subsequently 
by Launay^'^ . The main assumptions and apo roximations involved are 
the same as discussed in chapter II and the Cartesian Co-ordinate 
system is identical to that described in the last chapter for the 
body centred cubic lattice . I general substitutional injnirity is put 
at the origin whose mass is M + ^ M and which is coupled to its 
nearest neighbours through force constants +A V-j_ and ^ 

4-AV 2 . The total number of atoms is taken to be M. Vl anf? V 2 
are supposed to be unequal because if they are equal one gets 
results identical to those obtained in the last chapter. This is due 
to the fact that when ■^1 = , then for example the Cartesian 

X - component of the displacement of an atomr is not coupled to 
the Y or Z components of the displacements of the same atom or 
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Its nearest neighbcurs . Thus to make a transition to the scalar model 

from this, one has to set ^ 2 * 

The time independent equations of motion can be written in 
47 

the matrix form ^ ( ve shall write " In this chapter) 

DU - M C-l 2 IJ = 0 ( 5 .1 ) 

A 

where D is 3N x 3N matrix and U is 3N - dimensional vector- 
shall follow the following partitioning. U is partitioned into three 
F - dimensional vectors ^ oC = 1>2,3 such that Uj_ has the 

X-displa cements of all the N atoms as its elements. Also D is 

A 

Partitioned into the corresoonding N x N blocks each of which 

is s cyclic matrix. Using £ , m, n defined through equation (4,1a), 


i.e 


R = ( £ + -g—) ® £ 1 ”g“ ® £ 2 ® £ 3 

Vl -Y2 \ 


and defining the two constants 

A 

the elements of can be written as follows, 

D^^ { irm , ) 


^ ^ L\i ^m',m 8n' ,n 6'n' ,n|^Sm' ,mtl 

S m» ,ffl ^ ij ’-Sn,n»^Sg^' , ^4 1 ,m , 1 

* 8 V » 2^- 1 ^ m' ,m + ” + 1 6*^n’ ,B-1 

S m' ,m - 1 , 2+ 1 5m» ,m - 1 ~Sn' ,n+l 

^2' , £-1 S m’ ,m + 1 S n* ,n - ^ 


( 5.2 ) 


for oC 


1,2 and 3. 
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Dj _2 ( ^ £ ’ni*n’ ) = Dp^ ( £ ron, £'n'n' ) 




I'r 

n»n’ ) 

/ 

1 

c 

,n 

5 n' m 

■ n’n 

r 

1 

1 ^ V 

ti 


21 


Ij "'' 0^1 <Sn'n-l 


f 1 O 2i' m - 1 



(5.3 ) 


+ X r 

J2.t f-1 O m‘m + 1 6 n'n - 1 

^ j' J 

( £nin,£'in'n' ) = ( Irm, ) 

r ^ r 

2 ^ nj n m'^Q+l ^ £’ ni'm+1 Sn'^n-1 

S' ^ n'^n'<*l^ (^£'^£+1 (5 n^n+1 

^ ^ L' L~^ Smjia+l 5 n'^n-1 ^ ~ S'n^n Sm' 




and 


^25 ^ ^ £'ni'n' ) 


m^m+l 


m'm-1 
(5.4 ) 


( £mn, X'm’n’ ) 


\ ^ ^ 

1^2 S ^ n^n^<S m’^m+l 5'^* ^njn-l 

'*’ S m'^m-l S m'^n+1^ £^+1 S m^m-l ^ n’^n+1 

^ S'j^’ £—1 ^ m^m+l c)n'^n— 1^ S'n^m— 1 

* ^Jt’ £-1 S^m’m +1 ^ 


(5.5 ) 


The solution of equation (5.1 ) is denoted hj U© * 
%bstituting i 

(r) = exp (i k - . R ) 
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in equation (5.1) leads to the eqaations 

- MUD 2 =0 

. 1') 8xp/i k.CEi - R 


( 5.6 ) 


iviie re 


V R' ^ 

(defines the olenients of the djmamical iratriy (J* • This equation 
(5*6 ) becomes ^ 

"v = M V (5*7 ) 


vnich is an eigenvalue equation. This gives three eigenvalues 

p 

MC* 3 ''(k ),3 = 1,2,3 which are the three frequency branches and 

j 

the corresponding eigenvectors v denoted by e^^'^(k ) with 

components e^P (k ), =1,2 and 3 ; are the polarization vectors 

cC 

for phonon propagating with wave vector k . For the model under 
consideration, the (j" matrix has the explicit form. 


cr = 


8^^ ( 1- ) 8 jSg SPg ^^^2^3 




^ 2 ^1^2^3 


with S-’ = Fin 


h 




B/3p 8|3^(1-C^C^C3)^ 

kn- 


(5.8 ) 


and ^ = cos 


The eigenvalues K (a. 7 (k ) are found from the roots of the cubic 
equation, 

X 1 t 0 (5.9) 


with the roots J\^. . =1,2,3 related to the eigenvalues 

^ *3 


as 




(5.10) 


'(k) are 


Then the eigenvectors 
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\C 


A A 


. . s s + s s c 

L '2^123 123 

' 1.^ 


3 ^ ^1^3 - =lV3) - 8^, 0^S^S3j/D. 

- Y/s) -AjIa - \%^iA 


xA (as * _ 

(5.11) 

v/ith j — 1,2,3 and D. am proper noniiali 2 Stion factors. 

t! 

The polarization vectors satisfy the orthonomality and closure 
lelatiois , 

( 5.12) 


> -S.f 


end / ) 

J oC - 


( 5.13) 


( i ) ^ f 1 'i 

and products like e^ (k) e ^ (iS ) transiPorm as k^ k^ under the 

syminetiy operations of the point group of the crystal xmder which 
2 

LO . (k ) remain invariant. 

'’'’'ow let the substitutional impurity be put at the origin of the 
co-ordinate system. The equation of motion will now read as. 


DU -MLO^U = PU 

whose general solution will be 

A A 

U = Uq - ^PU 

where the Green function G ha:S‘ elements ■ given by^'^ 


(5.14) 


(5.15) 


CJ 2 _ (.)2 V 


**“ 3 


2 . 

In Particular^ ^ ) is diagonal in c6, ^ for 
the present model due to the transformation properties of 
e^'\k) e ^^\k) under the ^sanmetry operaticns of 0^, the point 
group of the crystal under ccftside ration. 


( 5.16) 
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Thus 
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, 0,0 ) = 



1 

2 ""TT 


^ -IjJ . (k) + 

tl 


(5 .17) 


which shows the equality of the three diagonal elements as well, 
^ 5.5 The Isotope T)efect : 

From eqaation (5.14) one gets 


^ /Sj3iR}R" 

(5.18) 


Now let the substitutional impurity be an isotope of the normal 
atom. Then Ar, = AVi = 0, and the perturbation matrix P has 

non - vanishing elements corresponding to the defect site idiich is 

A A 

taken to be the origin. Noting that G and ? have been partitioned in 

A 

the same way as D one has ^ 

(£',2" ) =AMl.3^S/3/i’ (6.19) 

Substituting this in equation (5.18) one gets 

V® ® ’ -’2^ ''/S ^2) (6.20) 

A 

Putting R = 0 in this ^ one gets , 


n , ( 0 ) 1 -- —5 . (0 ) 


1 + G (U) ?0,0 ) 

2 

where the fact that G ( CJ ,0 ,0 ) is diagonal in oC and /S has 
been used. Substituting equation (5.21) back in equation (5.20) yields 
the result . 


(5.21) 


= V® 


A M OO 

j 1 + AMO)^ G( tO?0,0) 


G^^(00 ,R,0) (£) 


(5.22) 
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As discussed in Appendix ^R,0) can be expanded 

for large R as 

" "Z ^ ^ exp (i I^. R ) (5.23) 


Inhere each term corresponds to a frequency branch denoted by the 
index 3 . Substituting this eouatibn in equation (5.22) one gets 


V ' V'®’ "’2 ' 4(^4) (ifep- B 


) 


I'P R 


where 


AMtO 

1 + AMt0 2G(c0f 0 , 0 ) 


1^ 


(5.24) 

(5.25) 


lere the index i refers to the polarization of the incident phonon, 
'^hus we have the solution for the scattering problem in (5.24) for 
the isotope case. One can fix the incident wave by choosing TJ^ 
suitably. It is clear that for a given incident phonon with a given 
frequency and polarization, the scattered wave consists of in general 
three polarizations. Thus an incident vibration propagating in a 
certain direction with a given frequency and polarization is scattered 
from the substitutional impurity and at a large distance R the 
resultant vibrations can be analized into three distinct waves 
travelling in general in three different directions with three distinct 
dispersions corresponding to the given frequency of the incident 
vibration i,e with throe different group velocities. 


k 


, 00 





o . ui 


Incident 


Impurity 
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We shall call the three final polarizations as the three 
polarization channels. Equation (5*24) permits one to separate the 
Contributions to different polarization channels. Thus , if one is 
interested in the scattering to the firf^t channel ( with j = 1 ) only^ 
then this process is described by 


(S.) + >ko ) exp (ik^.H ) (5.26) 

~ Kfi n 

Similarly for the channel j = 2 , one gets 

^ (k,k^ ) (5.27) 

NiR 

and so on. Then these scatteilng amplitudes can be used to find 
contribution to the total scattering cross section and mean free path from 
phonon scattering into different channels. 


However one can employ the optical theorem to find the total 
cross section directly. This involves the imeginaiy part of the 

forward scattering ■■ arplittjde i which for the present case becomes. 


d (k,k) = 

3 


■ V , 

> g (k ) 

4sMC0"^G( bO fo,g)Z__ P 


e^'^'^Ck ) 


where k is the incident wave vector and j is the incident 
polarization index. 

In particular , let the phonons incident along (lOO) 
direction be considered. Then k = (k, 0,0 ) . Further let them 

be polarized, say longitudinally with frequency given by ^(k) • 
Then we have 


(5.28) 




( 5 . 29 ) 
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and 




Am lO 


1 + AmoO^Cu) ^,0,0) 


(i£ ) 


( 5.30) 


In the limit of long waves one gets for this incidence, 

«?"■ ^ -L'S '«.*>] 


4Tr 

~k”~ 


AM lO ^ g^(k) 


U' 


(l+ AML0^G( tA ^0,0) I J 


2TT‘ 

But g (k) = — 

1 /3,v- 


16 TT^ 

M V*U)^ 


AMlO^ Im G(ua ?o,n ) 


(5.31) 

(5.32) 


where 00 


^l/M 


Also Im ,0,0)= + 


Im 



3 »M . i e 

k 3 ^ 


Im 


3 MV 



TT to \ 

6 M V* ^ 


k 


60 ^^ijO^ik) + ±e 
3 3 


( 5.33) 


for small oo . 


(5.34) 



X cos ( ^ 


(5.36) 
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~ (27)2 gin 0 cos^© co^(n sir? (p 

where cos = L 

(cos^© + sin^0 C0S^(^ sin^ (j} 

We have also to use the fact that for , 

k = (k,0,0) , = orCO^ = 


2 ^2 
^ Ca) o 

' i 


(5.37) 


so that 


k = ^ 

‘ ’ "~u) ■ 


Thus finally we get 



_4 , A M .2 / (>3 ^ ^ 

IT ^ M ^ ,0)| 2 


in units of a^. 


where c^^ = 


"4j2Tr® - 

. 

b = 

^__£r_ 1 

_ 6V*8 _ 

1 

_ 24 X 16 “FT J 


b. 


(5.39) 


The same result is obtained for all initial polarizations with 
(loo) , (OlO) and (OOl) incidences. This result goes over to that for 
the scalar model given by equation (4.39) of chanter IV ( one has to 
set AV= 0 in this equatioi to get the result for isotooe defect 
case) if we set /2 2 =0 i,e = y 2 so that ^ ^ ~ ISTfx 16 \JJi 

from equations (5.35), (5.36) and c^ is unity. This transition from 
the present model to the scalar model by adjusting the force constants 
suitably is interesting to note. 

The expression for ^ ^ above shows 'Rayleigh type of 
scattering with possibilities of resonance at xdiich CO~2 ^ 

The remarkable point that emerges from this is that although the 
dispersion used here does not go over to Deign* model dispersion in 
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the limit of long xraves like the scalar model, wc still get the same 
result except for a numerical factor. Thus our lesult diffi^rs from that 
obtained -by using Debye model by a numerical factor. However, the 
difference lies in the fact that this ntimerical factor depends upon the 
direction of incidence the polarization. For example , for (llO) 

direction of inci>'.ence , c^^ in equation (5.39) must be replaced by 

^2 “ (^2 -/ typical initial 

choice. This reflects the deoa .ture of the present model from the 
Debye model in the limit of longwave s . 


^ 5.4 General Defect ; 

Wow we consider a substitutional impurity with changes in mass 
and force constants, let M + A M be the mass of the impurity atom 

which is coupled to its nearest neighbours through force constants 

A 

')f + and 'Y ^ + A iT n • pertuibatim matrix P can 

J. JL 

A A 

be partitioned in the same way as G and D as^ 


r— A 

A 

A 

P, 

P. 

P 

11 

12 

13 

A 

A 

A 

P 

P 

P 

21 

22 

23 

A 

A 

A 

L 31 

^32 

P 33 


and each of the N x N matrices have non-vanishing 

partitions p^^^ in the r - space j the space of affected lattice 
sites. Then one has 

A A A A .._to7 

Pll = ^22 = P33 ^ij " ^ 

with p^^ (0,0) = .Am<a^^-8A|3>^ , P^^ (0 


~ PllCSn^S) “ Pll >Sn) “Afi > 

Pij (2 /Q. ) “ Pji »Sn ) “ 0 > i j 


(5.40) 
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Pij (0 > Sn ) = Pij (Sn ) = ± ^ j3 2 ’ i j 


where the sign depends on i and j and the nearest neighbour 
whose lattice site is given by . 


Here 


+ 24^2 ) 

3 


and A 2 


( 4^1 - 41 ^ 2 ) 
3 


(5.41) 


We follow the partitioning already described and define 
A A 

S=PU,So=PU^ "to scattering equation as 


A A 

S = ^ - P G S 


whence _s^ ^OoC. 


A 




where is the only non-vanishing part of 

in the space. 

If vre define a 27 - dimensional vector , 


and 


(5.42) 

(5.43) 


lies 




3,1 


A 

s. ~ 


-2 

f 

P = 



1 s.. 










A 

A 

A 

A 


^11 

%2 


and g = 


A 

A 

A 



"21 

"22 

"23 



A 

A 

A 



%i 

®32 

^33 

then we get 








A A 


s 


io - 

P g S 


whose solution is 






A 

A 4 ._1 


s = 

(I + 

Pg ) 

®0 


A 

A 

A 

^1 

^12 

^13 

A 

A 

A 

^21 

^22 

^23 

A 

A 

A 

P3I 

P32 

^33 


•(5.44 ) 


(5,45 ) 
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^4- A A 

= 7 ' M V 


( 5.46 ) 

where the V - matrix found group theoretically block diagonolizes 

A A A ^ A 

A 


AA A 

(I + pg )~^ to M. 


% canvas in chapter II, write M as ll = ^ corresponding 

to different symmetiy blocks and get 

A A 




_s - V M V ^ £ r 


(5.47 


so that 


From these we obtain expressions for _s^ whence S ^ are 
known and then 


gives 


2 = Ho - G s 


(S) ~ S^Sn ) S^(£n) (5.49) 

For large R one can expand ^ , Sn ) asymptotically 

to get, 

V(S) ' (io > f (si sj ) exp (U^-B) 

j I R f 4^ 

where k and i specify the wave vector and polarization of the 
incident phonons and F (k^, ^ ) is some trigonometric function of 

k and ^ . Then one can use the optical theorem for the total 
scattering cross - section. Ijt addition to the asymptotic expansion 
of Green function , one needs the elements g ^^ ( ^>Rn#Rn ) 

and the V - matrix for explicit evaluation of the scattering 
amplitudes F (k* f |[^ ) and their partial wave expansion' . 
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proceed to discuss these in the following. 

To find the T-.^trix mat firt the rep„aentation 

! of 0^ io the aoace of the displacemente of the impurity and its 
nearest neighbours. Por the .„del under consideration , this space 
13 tvsrntyseven dimensional. Howeyer it is the direct product of a three 
dimensional space consistir* of the Cartesian oo^onents of a displacement 
and a nrne dimensional space consisting of the position co-ordinates 
of the affected lattice sites. Hence T is a direct product of f j , 
the representation of 0, in tee ( x , y , a ) space and P ^ „hi=^ ’ 

IS the sams as that for the corresponding scalar model rfrioh for the 

present case is P =21 +p ^ ^ 

2 Ig ^lu ^2g •^2u * we get the 

result also obtained elsewhere ^ ±n a different way, 

T' r, ®n 

"I'l ® "'ig * ""lu + I-Sg * -f-ar ) 

“ hg * ’ ^lu * ®g * ^ ^2g * ^2u ''an + % t (5.51) 


where we have used the results , 


The above equation ininediately shows that in the partial 
wave expansion, we will have S, ?, ^ r, 

of waves where the noaeftclature is that depicted in Appendix XII. 


^lu(X>\g 

= hu 


^lu® hu 

= hg ^ \ " Ejg * F^g 

(5 .52) 

i 

"lu® F,. 

= ^2u " Eu * 

1 

’'lu®J>ar 

1) 

( 
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I 

.9 


f 


( 


The sum of ^2 '’ 'wave and that of and 

£"2 is the P wave. We have distinguished between and II 2 because 
they arise due to different irreducible representations of 0^ . 

Av 

The Y - matrix can be evaluated by using the projection 
technique desl^ribed in Appendix III with the aid of equation (5.51) 
for P , An alternative method will be to use the Ciebsch 


Gordci Co - efficients for the direct products given by equations 


(5.52 ) and use the known basis functions for P' 

A 

columns of matrix Vg given in Appendix III along 


the basis functions for 



simply 



which form the 

2 

with the fact that 



This method is elegant and useful. The resulting matrix is depicted 
as Vp in Appendix III , the columns of which are eigenvectors of the 
symmetry types in the following order , 



P , S ^ Dg / Pp > Pg > ^ and H. 

The rows are labelled in the same way as those for the 

A A A , 

matrix p . In fact Vp can be partitioned in the same way as p to 

give , 



and a conparison with Vp depicted in Appendix III gives ihe 

A 

structure of these 9x9 submatrices V^^ . 

For example , 
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V. 


11 


with 


1 

0 

0 

0 

0 

0 

0 

0 

0 




0 

0 

0 

0 

% 

0 

0 

0 

C 

0 

0 

0 

d 
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0 

d 

c 

0 

0 

0 


0 

0 

-d 

c 

0 
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(5.55) 


oimilfir ly on© can write down other subrnat rices. 

It has already been pointed out in chapter II that the 

stiucturo of the block diagonal form of p , g , (I + pg ) and 
, ^ 1 

( I + pg ; con bo known from the structure of the reducible represen- 
tation P , Thus from the expression for P given by equation 
(5.51 ) and tho order of columns of the matrix 7p described, we 
expect tho follovdng block structure given in order. 


ThrfJo equal 3x3 blocks for P— wave. 

Ono .single element for S - wave . 

Two equal single elements for - wave. 

Thiro equal 2x2 blocks for Dg - wave. 

Ono single element for - wave. 

Throe equal single elements for F - wave. 

Throe equal single elements for G - wave , 

Two equal single elements for H - vrave. 

2 

The structure of Rji , ^ ) can be written 

down using equation (5.15) and noting that ©'■^ (JE ) 

transform like under inversion and group operations. 

Also wo note that (k) (k) is real and If we set 

o( P 
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S'. 


, Si ) = (a, , a ) 


and 


Sn-Sn-M^g^pCCO^, Rn,B^) = \^\ ) 

This leads to the concluGion, (we suppiBSS tO ^ in tie following): 

(5.56 ) 

’ Si ) = ^ - ft. ) (5-5’ ) 

T’or the model under consideration, 1 can have five values 

which appear as subscripts 0 to 4 in the following which in terms 
of tho components read ( C , 0 ,0) , i ±-7 >±7,1^^) > ( t 0>0) 
and its cyclic permutations, (il, + 1^0) and its cyclic permutations 
and ( +1, +1, + 1 ). These are the same as those for the scalar 
model. Inserting these values of ^ into the expressions for 
) we get the following different Green function 

elements, 

^dd ^Sn ,Sn ) > o( = 1, 2, 3 we have 


^Do = 


N M 


■zz 


I 1 


^ (k ) - 

2 ~ 


3 N M 


I 




^ (jO^ - 6L)^(k)-i£ 


( 5.58 ) 




N M 



e^^(k) cos cos cos -|k. 


oC '^•4 
2 TT2 


5 N H 


U> ^ - 00 . (k ) - i 6 


cos .p-kj^ cos 7 ^ cos akj 


(jL)^>U)^(k) - i€ 
2 


(5 .59) 
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"D4 


^ ST* 'V ^ ^ ^ ^2 ^3 

NM Z^Zw _ 4^2 ^ ^ ^ 

£ 3 2 ~ 


hr 

Jf<'>s] I 

•# 


g. 


1 cos k^ cos kg cos k^ 

T¥T L L ']Jz7Zjz~iyZTi 

i i j “ 

^ TV® V 


D21 F M 

k j 


u) 


u)2 (k) - i e 

2 


oC = 1,2,3 


g 


D22 N M 


^(i) (i) 

1 “S;^ r~ (ii) (k) cos k 

•^4. Z* (jt)^ _ (k) - i 6 

k 3 3 


g^ 


D31 


■wIE 


k J 


Ji, P> = 1,2,3 , 

*( a) ( i ) 

e^^'^(k) e^ (k) cos k cos k 

’"to 2l’'J>2”k) - i e 

j 


and 


g 


D32 


for c/. ) ^ ^ 1>2,3, ^1 

Zj Zj 


( 5*60) 


(5.61 ) 


(5.62 ) 


(5.63) 


NM Lt Li (O 2_ ^3 2 ^ ^ 

k J j 

ot , ^ 

such that ot / /3 , ^ 

arkl 


(5.64) 


4?f' ' , 

4i^ 


* • ^5 • « 



d 



hese aie the independent elements cccurLng in ) , 

^22 (Sn , R^i ) and gjg (R^ , ^ ) . 

/ o(^^, </, ^ = 1,2,3, we get 

a - 1 ei^'^(k) e 1 -r^ ±^^ 3 ) 

%1 - -^JT" y > 

k 1 ^ ^ - O)^ (k) - i € 

(5.65 ) 

= ^ V T ® ft ^(is) ^1 i ^2 ± ^3 ) 

nm AZ. "^t'TTp ( 5.66) 

k. 0 ^ - u) 2(k) _ i £ 

«i 

"« A Z ~-t;)T:-;3-r--:7- ca.s7 ) 


Iho subscripts D and ¥ lefor to diagonal and non-diagonal 
Partitions of g. The first index refers to the j | values and the 
second index refers to the different types arising from the same 
llkl value . 

With, these one can block diagonalize g and p and hence 
A 

( r + Pg ) whose inverse can be found and then s ^ vectors can 

— 0^ 

bo calculated . 

The resonance denominators are given by, 


F. 


v 

^ ^ ” %4 ^ ” (^D21 “ ^D3l ^ 2( gjj22 - gj^^) 


- 2 (gK3 + gN4) J 
A rs 


Fp = det I (I + 4p ZSg ) 


(5.68 ) 

(5.69 ) 






• « B6 « « 


where 


A 

A» = 


aMU>2-8A^^ 8 cA|3j 16^p2<i 


and 


A 

A 


3 


16A^jd 


^DO 


Api 

0 


~ A ^ 


8 c 


®D1 


16 d g 


m 


1 %0 ^ ^ ^DPI + 


M ^ Sd 21 "■ ®D31 Nr2(gN3 •*■ %4 ) 
+ 2(gj)22 + gD32^' J 

nT^C gN3 ^ gN4) l[gDO ^04 “ gD21“ 

“ SN3 %4 ~f 

•iJ 

whom c and d are defined in equation (5.55 ). 


^Dl " ^ 


" ^D4) - (^BZl - Sd3i^ ^^^D22 " %2 ^ 


D2 


Ai-api 


* (g]V3 gu4 ) J (5.70) 


1-Ap^ 


2( 4 (^l) ( % - gH4 ^ 


^P1 ~ 1 


F2 


^^D21 ■ 

■ %3i> - '^m * 

%4^ 

II 

^ ~ ^®D21 “ ^D31^ ” 

D22 ~ 


%4 ^ 



( 5.71 ) 

^%21 ' 

^ gD3l) ^(SD22 ®D32^ 

” ^^^N3 

~ %4 


(5.72 ) 

^%21 

gD3l) •*■ (gfis “ 


r , , 

(5.73) 
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04 , 



♦ ♦ 3*7 » • 


F(, = 1 - (Ajij,) 


(SdO ■ *D4' * ^®1I3 


(5.74) 


and. 


T^h = 1 


(Ap,j^)^gl(j + gD4' * ^«D21 ^ ®D31^ ■ ^ 


Sd32^ 


+ (%3 " %4 ^ 


(5.75 ) 


An alternetive way of finriing tho y yectors is to use 
tno 01 ol^.h - Gorton Co-afficients ani to Ite use of the rtsults for 
the scalar moael. This is easy to io with the single factor blocks 
For bigtor blocks the problem of indexing the Orson rtnction bakes the 

D ro CO c 3 c omoli ca t/ed • 

once the 3 ^ vectors are cs,lculatea , cno can firt 
F (ki , kl ) aopea’rihg in equation (6.50) . The expression for 

F (k^ k^) contains contributions from S , p, I>i > “s’ 1’ 2 

” -1 Hovjpver in the limit of long waves, 

0 and H type of waves in general. Ho . 

the F , F_ and H type of waves do net contribute . This can be 

^ ^ aor the scalar model whera in the 

snnticipatod by examining the resul 

nnl v the waves corresponding to 2 lu 

limit of long waVGS only the wa 

. Hence in the present case one expects to find contnbut 
contribute . Hence m me hx« 

from wavon corresponding to F^y^ -^Ig 

(5.76 ) 


= 2 F,y - - Eg ^ F^g - Fgg 


i 3 ^ 

„ +>(0+ contribute to F(k , k ) 

X Vh.«b<.-i aelv those waves that conux 
pvnd this loads to pree Y 

c TVma in this limit , 
in the limit of 1°^ , „ 

F 


u limit of long waves, xnus 

■ (S)/ i 1 ^ \ + F^^^(k^,kh ■*■ ^ ^ 

F (kSki ) = f 






(5.77 ) 
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where k and i specify the wave vector and polarization of the 
incident phonons and and j specify those for the scattered 

phonons when lie along symmetry directions. Using the long - 
wavelength limit expression for incidence along (lOCj direction, 

g a (kj) = —€—.Z2.L. = 

f T ^ 'TTr . , . K O 


TTxsA TrMU)2 

* n 


with (jQ 




the expressions for ^(k^, k^' ) 


M 


found to bo 


(5.78 ) 


(kS kj ) 


■(lTr.“)(-|;j[Z *<''■**> 

/ - 0 I’ t-£/ ^ 

£' 


(5.79) 


'' o P / 


with 


^2 ~ 


IZi 

= AmIa) 1^22^33 "A 23 ^ 32 ] 

^^^12^25 “AggA^JJ 

vAl^ 


(5.80) 


(5.81) 




^ FAei a 13 


A 


A A 

where /\= 1 +A_ , 

^ p g > 


A 1— 

* A 

the 3x3 matrices and /\ 


(5.82) 


g 


Sie defined after equation (5.69 ). 
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4 ’(k) J 


- k, 


o2 ®2 ' ®l(i£) ■*■ ^ ®^(l£)" 2 k 2 e 3 '(k^ 


2^03 






X 


X ei*'^’^(k^) + k^2 ®2^^'^(^) - 

(^ ^ \ f A g_l \ 

s STTFpj, J \ jB 1 / 

|^>'2 ®i(» + ki iitk)| /k „2 
ka e^'ca+kj jita ^ fk^ e;(j)(kj) + k„i ej'5>{ki ) 


(5.83) 


+ koi 62 ^^ ^(ki ) 


] ka e.^k) + fc. / k 


:o3®2^^'^(14^ ^ ^02 ) (5.84) 


c 


J u 



whe IB 


X 4 (^Do ^ ^ %21 *■ %3l^ ^^®D22 ~ ^D32 ^ 


and 




kj) =/'- — 

\ 2 7 


) / ^Pl 
TTfJ V 



^esCk) - kjGgHk) ^ X 


(5.85) 


^o? ® 3 ^^^(^o) - ^o 3 ® 


{j^os 


2*^’(kj )j' + /kj ijtk) - kj_ ijik) 


) - k 

o '' ^ol ®3 


y l^k^T ep^^'^(k^) - kn? e/^'^(kj' 


^ol ^2 ” ^o 2 ®1 ^ (feo^ 


(5.86 ) 


In all the equations above one must evaluate the quantities 
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^ ^Dl , , and in the long wave length limit . 

Ltt the /'incideny phonons be incident along the (lOO) direction. 


The n k. 


1 k, k2 ~ 0 ^ = 0 


anc* the total scattering cross section 


is given by the optical th 


eoreci* 


Im F(k^, k^) 


Polarization Here • 

I'ow F (k' , k ) with k and i specified as above leads to the 
following partial wave contributions, 


(5.87) 


^ , kb = - 


3 TT F. 


sttf. 


Tf Mu)% 

O ^ 


(h^+ 8h^) 


(5.88) 


(ii) Transverse polarisation Here e^^(k) = d^?. otSb-Z 

vis) = 0 = F (d ) 


.(f*) = 


TV Fp 


(h3_ +8h2 ) 


F^^2^ = ( 1. 


(fi) 

F^ ^ = ( 


2 IT Ft 


s-TTFg 


) (^hl — ) ( l-ASi Ap Ic' 
Pl 


) ( 

di 


(5.89) 


Mow we have for the Green functions in the limit of long waves 

Xg kg Xg bg 

Ijf gpo = ^2 , Im + "jg~ > I® §^4 = *£-+- “jg 


(5.90) 


1 . 


”1. 
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with 


z_. 


sin Q 6 9 ^ 

f (e /f) t ^4 


( 5.91 ) 


'hcrp (0 ) is given by equation (5.35). 


( + Xgbg ) 


I® ^DPl ( ^2^3 

) , 

%31 

%22 ^2^^ 


®D32 

Im %T1 “ ^2 4 ’ 

^N2 ^ 

4x2 4 

txJTT 

where x, = - 

\ 

V - 


( 5-92) 


TTiJ 


6 M V U}^ 
0 


M Y"*' (x) ^ 


and b^ , i = 3 to 8^ are cumbersome integrals over angles like 

and bg described in equations (5.35) and (5.91). These yield 
mimerical values after integration with given ratios of ( — -- ) 

/3 1 


Using all these results and the optical theorem along with 


the fact that for (lOO ) direction of incidence, k = 


.2 _ 8 lO' 


one 


(A)‘ 


Uots the results 

'(100) 


O', 


(100) u /(O' 

2 


/A^I 



ITT/ iw *itf 


hp! 


m { ■. • (■ 


=4* “5 






-KrKlf 


^2) 

(5.93) 


and 


^ (lOO)T /j£ 
^ T ~ I .V 



UO 
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2 



(AB>2\ 

I 1^2/ 


? 


( 5 . 94 ") 


Here T and L refer to %e transverse and longitradinal polarizations 
roapoetivclv and the numerical cfuarrtities a^, i-y= 1> ...» 10 are 
appropriate combinations of b^, i = 1, ••• > 8 and other numerico i 

factor?' . 

Pimilnr erpiessxons can be derived for other inricont 
diTT^cxions and polarisations. However the partial waves contributing 

for different incident directions and polarizations are different . 

Table IT bolow gives the result in the limit of long waves. 

TABLE II 


Pi TCoti on. of 
Incidcnco . 

Polariaation 

Normalized 

Polarization 

vectors 

(transpose) 

Partial wav 
es contrib- 
uting to 

CT T 

L 

Partial wa'^' 
es interfe- 
lingin the 
expression 
for /C'(t^^ 

Bemarks 

lOO 

Lengitudinal ( 

Transverse 

Transverse 

;i 0 0 ) 

( 010 ) 

( 001 ) 

S,P,Di 

P 

P ,D2,G 

S-P,P-D^ 

P — Pg 

P-P2,P-& 

The two 

transverse 
branche s 
are degenerate 

110 ' 

[jongitudinal 

Transverse 

Transverse 

1 0) 

(0 0 1) 
JL-(i-i 0 ) 

v'2 

-Ml 1 1) 
M(i-i 0 ) 

JM 1 - 2 ) 
fe . 

S,?,D2 

P,Di,G 1 

S-P 

p-D2,p-g ; 

?-D3_,P-G 

Transverse 
branches are 
non-de ge ne rate 

111 

Longitudinal 

Tran.sverse 

Tranave-rse 

S,P,I>2 

P,5i,D2,G 

S.P,P-P2 

P-G 

P-Di,P-D2, 

P-G 

Transverse 
branches are 
degenerate 

.i J 
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It should he noted that even in the limit of long waves , the 
number of partial waves contributing to the total cross section depends 
upon the direction of incidence. This feature xs a consequence 

significant departure of long wave behavior in the present model from 

^ nf pd to the 'behavior 

the usual Debye approximation. It is in fact rela 

of the dispersion which in the limit of small k becomes. 


.(k) = 

3 


M 


Its 

2 M 


'k> 


with 


for '3 = 1^2,3 
2 


27 ( k^ k| I?! 

^ 

( ( 1-2 1.2 , I 


( 5.95) 


( 5.96 ) 


(kf + if 4 * '"s 


This shows that the .onstsht f«que..y shrfe.es in this 

long v^veUngth lisat ate act sphere, as wouW he the ease if Deh^ 
^.el is ass«i. - nthough this aeperture ^owei itself wp for the 

oeee of ieotoos iefsot in the ».rfcei factor as 

X nf « general substitutional imourity quite Slgni i. 

affects the case of a genera 

rtvnp oath will contain interference term 
Por ox, mplo the inrorse man free path 

. aifforant Partial wave oontrfhutions aepenfling upon 

hotwnon different p, 

X- shf -ineidfince and the polarization o 
direct Lon of m - interference 

ThU is Shorn, in the ahovetahle.lt is to heated that an 

oocors hotween waves of different 'pamUes- . 

-Lv. rross- section one finds 

TTrora the expressions for e 

f mass affe-ts o^y the P -- oontri^tion Pereas 
that the change o ^ the waves. The low energy 

tho change ift constants af ec s ji^^yleigh type 

SI. henoo long 




Oi scattering in absence of resonances. However the presence of 
resonance denominators include the possibility of resonance scattering. 
A comparison of the results obtained here with those of the scalar 
model discussed in the last chapter shows that there is a change in 
'parity' of the waves in the two models- For exanple the, S — wave 
of' th"^ scalar model goes over to the P - wave in this case and so on. 
"Phis result is expected by simple group theoretical arguments. 

* 5.5 Nummary i 

Tn this chapter the model considered may be claimed to be a 
atop forward in the approximation to the realistic case, compared to 
those considered in the last two chapters. As expected one gets some 
results that are qualitatively the same as those for the scalar model. 
Among these are the modified Rayleigh type of scattering at low 
frequencies with resonance possibility and the interference between waves 
of different' parities ' in the expression for /V . One of the 

differences lies in the consequences related to the- validity of Debye 
r'Jproximati.on in the long wavelength limit and the shape of the constant 
Tt 'qiu roy surfaces in this limit. It turns cut that in this limit 
1 ho constant frequency surfaces are not spheres in general and hence 
I'.ibyo ap p roximotion does not hold. However it is remarkable that for 
+>0 irotopr cn£3c this leads to the same result as that obtained by 
using Dobyu approximation except for a numerical factor which depends 
upon tho direction of incidence. For the case of a general 
r.ubstltuti™l impurity it has the effect of making the long wavelength 
total cross - section deoendent upon the incident diiection unlike 
tho scaler model ^eie it becomes irdependent of incident direction. 
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”0 h,v„ also aiscussoa 00. can go ..er to the ecelar n_^el from 

the model under consideration by adjusting the force constants suitably. 

the other non feature here is of course the dependence of tie 
ecatternng process on the polarisation of incident phonons. It is 
rluuedy apparent that Introduction of polarLsetion leads to a groat deal 
. f rathematical complexity 1„ the problem. Tie soacterinc process 

quite inyolTCd physically also. Compared to the scalar model, 
tho cormsponding partial wares appear, as exeected,to have changed 

their 'parities' . Perhaps one can safely contend that these additional 
features mst be present in a realistic Case. 


CHAPTER Y1 


COIlClHSIOn 

• Sf 

he r.‘fiulto obtained in the present study can be susr^arized as 

illo'js , 

The scattering cmss section and the nmber of partial waves 
'rtr. butiiig to it aepend unon the direction of incidence of plionons, 

U'a,s <? hi biting the expected non-central nature of the 'potential' due to 
. sub;' bituLi onal impurity m a lattice (Chapter IV, Table I). 

Tnn resonances in tlie partial x^ave amplitudes and the behavior of 
i'UC ntlciiiig cross sections near resonances have been investigated in 

ti' . '^ho tcsuIijS confirm those obtained by other I'/orkers. The 

ring of" both the acoustic and optical phonons in a diatomic simple 
' nl'ic .1 ai<ticn has boon studied and expressions for resonance denominators 
obLfiinc'd explicitly. 

Tho expressions for mean free paths of phonons due to scattering 
I ‘oiu ! ul utiintioral impurities have been obtained for small concentration 
I'- i m, n iri t.ic-'S and phonons with long wavelengths. These expressions are 
iiuwn to contain intorfcrenco terms between i/aves of different 'panties'. 
T'loir behavior nnar resonances has been studied. Hie similarity of 
■ for Iho throe cubic lattices has been discussed. 

Tl\n offoct oL’ polarization of phonons on the scattering process 
h;if! t/itin studied by using a suitable model of body centred cubic lattice 
■n { h'lptnr V. It has boon pointed out that the results obtained for 
ibis c^usc can bo baleen over to the case of the so called scalar models 
by nn ai)propriate choice of force constants. The Debye model is shown 
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^ for the exact description of the scattering process i. 

f loni_, waves, iho departure from the Dob;-e model in this limit 

dependence of the scattering process on the direction of 
incidence (Chapter 1 ^ Table II). 


The entire study has been based upon the harmonic approximation 

and nearest neighbour interaction models. The effect of anliamonicity 

piionons in the impure crystal is rather complicated, and leads to 

50 

■; iinito life-time for each phonon. This woiold affect the phonon 
'loun free path^ and the transport coefficients which depend on it. The 
pn^sent .analysis is valid at very low temperatures when the anharaonic 
off 0 0 to arc small, and would be adequate to cover such phenomena as the 
Kinks in thermal conductivity curves obtained by many workers^^"^^"^^"^®, 


Thu other essential limitation of tho models investigated hero 

ho tlinb tho impurity is talcen as a substitutional one. There is enough 

51 

ox]aor iinrintal evidence now to indicate that quite often the impurity 

iibon (or molucule) goes into interstitial positions. The symmetry of the 
impurity 'potential' in such cases is surely more complicated than xjhat 
lias b'-’cn ;i;r!uniod in the present analysis. A proper theory of phonon 
scabtx:rjug from impurities must take such situations into account. 


A??E>]DT7 T 


GTOI FTOICTIONS FOR Sr AT, at? Mnn^T.Q 


The elements of the Green function matrix for a scalar model 
are of the form^ 


2 

G(a) , 


li, 




eocp j" ik. (R_r. ) "[ 

0) (k) - tO - 


(AI-1) 


2 

whore cO (k) is the appropriate dispersion, the components of ^ and 

Rj arc inoasiired in the imits of the corresponding lattice spacings. 

I hen k is diiao ns ionless wave vector and the k integration is over the 

first Brill ouin zone whose volume is V*. For the simple cubic lattice 

k integrations extend from - irto +T with V* = (2 TT)^, the unit of 

iiH5asurement being the lattice spacing a"'. 

For face centred cubic and body centred cubic lattices the 

S3anmel.r;^'- of the constant frequency surfaces can be exploited by doing the 

48 

fi tog ration throu{^ l/48 of the BriHouin zone , lAdoich in the face 
c<'>ntrod cubic case is the part of the zone enclosed in the trihedral 
jmgic defined 'by (001), (101) and (111) directions; and in the body 
centred cubic case is the part enclosed within the tetrahedrcai bounded 
by the pianos kg = 0, = Tf , and kj_ = k^ with their 

ranges defined T^y 0 ^ k^, k^ ^ and 0 ^k^^ TT • It is 
fdso possible to combine these wedges to form a cube in the reciprocal 
space whoso sides are of length iTT i.e, from -ETfio + 271" such 
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that integration over this cube gives four tines the result of 

gration over the Brillouinzone for the body centred cubic lattice 
cmd twice the result of integration over the Brillouin zone for the 
face centred cubic lattice. We choose this cube for integration 
because of the anclytic convenience it offers for body centred cubic 

and face centred cubic lattxces. llith this choice and a change of 
variable one obtains 


-hiT 


, Ri’ = 


-1 



exp J d\ 


"W 


(AI-2) 


whereas for the simple cubic lattice 

•+7r 

) = -=k 



exp I'ik.C^-Rp J d\ 




(Ar-3) 


[ho subscripts B and F refer to simple cubic ^ body centred cubic 

.'md face centred cubic lattices respectively. The asymptotic esspansion 

ir. obtained I’rom equation (AI-1) for large R = | R^-R. | , and this 

is discussed in Appendix II. Since G ( uJ , depends upon 

^ [I 

I Si“l- I } follows, we shall adopt the convention of writing 

tho Cartesian components of ^-Rj in the units of the lattice constant 

/ 2 

’a' of these cubic lattices. Thus for example G (c^) > >0) stands 

D 

2 

for R^-R. = ae^ + aOg in body centred cubic lattice, Gj,( a> A A A) 

stands for Ri-R. = a(e^+ i ®2 ^ centred cubic lattice and 
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so on. Hero are the unit vectors along the Cartesian axes. 

The Green functions w ixhoccur in the expressions for the 


sc 


Littering amplitudes are conveniently evaluated in terns of the 


parameters E„ = 3 E = 1 - ^ 

S ® 


oS 


Cj 

oE 


and E„ = 3 - 
F 


cO 


Cl) ^ 

oF 


so that the frequency bands are defined by the ranges -3^Eg;^ 3, 

-14^E ^1 - 1 ^ K ■ 3 respectively. These Eg,Eg 

and are dimensionless energy parameters. 

For the simple cubic lattice, the Green function elements 
oi'curing iii tiic expressions for the scattering amplitudes can be 
derived from the general expression 


G(p,q,r,Eg) = "X, 


-^T 3 

^ - d k cos pk^cos qlCgCos rk. 


L(2rrrJJ 

-ir 


E -cos k^-cos kp-- cos i £ 


(AI-4) 


A, ' - l/(M C0%) = - Vlar) . integral inside the 
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parf-inthesis con be written as , 


= i |^C(p,q,r,Eg) — i S(p,q,r,Eg) (A j.~5) 


wi th 

C(p,q,r,Eg) = 

and 


x >0 

r 


dx cos[EgXj Jp(x) J^(x) Zp^) 


J 

0 


L -OO 

S(p,q,r,Eg) = J sin(EgXj Jp(x) J^(x) 


(AI-6) 


(11-7) 


Ho, re Jp(x) is tte Bessel function of first kind of order p. The 
numorilal values of the functions C(p,q,r,Eg) and S(p,q,r,Eg) for the 
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values of fo n n^ /-i ^ ^ 

tabulated ^ ^28 32 ' 

tabulated extensively^® t„ +v ^ . ._^ 

limit of long waves one gets 

-ll-_ i03 


1“ I(p,q,r,E ) = 

Q *' 


nTp TTcO 


bS 3>r2 7r c03 

OS 


2 2 2, 


(P + 


q + r ) 


"here „o have enploj^d tte foUoaihg hsthcd. 

l(^) and FgCk) are two fhncticns of k, then 
uitegraL 


(AI-8) 


we can write the 


Im 


^ k F (k) 


^l(k) 

uning the formula 1 ^ _P 

X - i^ X 

Now W 2 can further write 


=irjd\ S' (ffls)) 

+ iJT C (x) . 


(AI-9) 






lv/i(k) 


(AI-10) 


nhere the Integration is over tto surface of constant k where k is the 
cojutfon at P^(k) = 0. All these are evaluated at the limit of long 
waves. Explicitly one has 


^’(O^OAEg) = 

S(1,0,0,E^) = 

<d 

sIFfFa)^ 

oS 

c> 2 

>r2'7ra) ^ 

oS 

3X2 TTu)^ 

oS 

C(1,1,0,EJ = 


2 A)^ 

S 


s/awb^ 

oS 

C(2,0,0,Eg) = 

--J±— , 

zfFfjs 



ZTf 

oS 


(Al-lla) 

(AI_llb) 


(Al-llo) 


(Al-lld) 
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For the body centred cubic 1 'itt-i + 1 , o 

attice, the Green functione used 

xn tno resonance denominators 


are 


%( XI , 0,0,0) = g , O 

oB 


4- IT 


L(2Tr )' 


d\ 


= 2 


V ^ = E g -0 

B oB '^B 


E^- cos cos cos k^_ i^ 

"bo 

“ 2 Ab J }%) [j„(2E^x) . i Yj2Ej^x)^ . 

^ J (AI-12) 


(AI-13) 


V ^ , 1,0,0) = g =;^ 


2B ' kg 


■fTT 

Yr 


d k 


cos 


2k^ 


^(^rrr UK 

--•ir ^B “ \ ^2 V ^6 

'l>o 


] 



(AI-14) 


%(0 ,1,1,0) = g_ 


'SB 


~r ® 


cos 


cos Bkg 


cos k^cos k cos k - i ^ 


] 


■&0 


r 2 


:yta J^(x) Jj(x)[j„(2E^x) . i r^(2E^x)j 
0 


- S^^foBR i GET 
^'^BL. 3 3 


(AI-15) 


Gb( CO* ,1,1,1) = =8 8 


g _ 3g - 3g 


B^oB B'^B oB *"23 3B 
Ab[°“> 4-^®\] (Al-16) 



. . 3 )3 . , 


wn.ro _ i/m l/Sy , J and Y are Bessel functions of 

i '^rt ,nd second kinri . n-.e quantities QZa and GBI. defined above have 
bo.n cuiputrd and bh 'ir nu:nriccl values are availatdo , lie distri- 
h'jbion function of the frequency (squared) for th^ body cen led cubic 
1 .- 1 t.\co, wliLch equals GBI^ encepb for a nuaericnl fac^or^ her a dogarithnic 
s ' iyjlo"ity in tli'^ trddlj ol the la' id i'"' tniu scalar --del. 

n tie 3ou frequen-y (loa ■ urvel^i ■ ii) lir.lt o. ctxi ioT< o.r the 
noohod -uidicatcd in equations (aI- 9) and (AI-10) to get 


Jn 


" B 1b’ ^ "iB = \ (Ub" Ib’’ I 


3 ir) ' ^ = 115 " XsCaiB - (Al-i?) 

.3. 


wb Ml a 




rji'^D ; re used in the derivation of the long wavelength lii3i+ of the 
S( o btci'.ing mpbitud'^s, 

'"imilarlj', fei’ Lhc free centred cubic latcice jre obtains 5 , 


-t-fl 


fi 

V Ct) ,0,0,0) 


pv ,u,u,u; -- -- Ap 




d"’c 


E^-(cos cos cos k^.os k^- 
cos k 005 k )- i-^ 

(AI- 18 ) 
(AI- 19 ) 


■r 






G ( O)*", 1,0,0) 


= ^ ^ J 


V —3 

^ M 27 r) 




d k cos 2k., 


'TT 


E -(cos k- cos k + cos k_cosk i 
F 2 ^ o 

+ cos k cos k ) 


(AI- 20 ) 
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-t-T 



4 1 

If 


-(2F)^^ 

i 


d ic 


cos cos li^cos 2k^ 


and 


d ( a3 

p 




+ T 




f{f 


U^T)' 


J. 


E -(cos k cos k + cos k cos k 
^ 2 2 2 

+ cos k^cos k^) _ 

(iI-21) 

2 

d k cos 2k. cos 2k 

— P 


_ ®p~(coo k^cos k + cos k cos k_^ 

’ B * 2 


-] 


+ cos k cos kj - i ^ 
(Al-22) 


where - - 1 /<,m - l/isy . The long warelength Unit for these 


nisy 


"of If' ^iF ^ ^ ggF " ^ ^F^°iF“^^2p^ 


Tn = 4 a.^), m g = 4 - 


“' 61 '' 


"F IF 2P^ 




(AI-23) 

Wi 1) T’O 

If = =>2, = aM k = 2 cO/ 

As nontionod above , tables with ntimerical values of the integrals are 
availalxlc. The nunorical computations of these from the foms given 
above happen to bo quite tedious and tine consuming. Becentl^ Ilahanty^ 
h<is suggested the use of the nsthod of Fourier series esqjansion of the 
Grt:!on flaaction. This simple technique is Jsased on expanding 
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^(fij x) 


sin(n7rx) 


n 


with a (R) = ^ 
n-' v" 


r 


J 


V’" 


d k exp(ik.R) sin(n7j- E(k)^ 


(AI-24) 

(iiI-25) 


wiierc E(k) M Ct) (k) and x is dinensioialess 


energy such 1diat the bard 

is always defined as -r ^ ^ / s 

The real part of G(R, x) is given by 

Hk) 

In G(R, X' ) dx» 

(AI-26) 


he G(R, x) = 

ir. 


X’ 


Ihe use of this method leads to considerable saving of 
computational, time. It may be noted that the quantities numerically 
computed do not involve mass or force constants so that these very integrals 
c'an also bo used for Groen functions encountered in the problem of 
scattering of spin waves and band electrons frem a substitutional inpuriiy. 
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iiPPEWDIX II 


iSlMPTOTlC EXPANSION OF GREEN FONCTION INTEGRALS 


A typical Greon function integral can be written as 

> £) = S G^^ ( (AII-1) 


where 


,R) 




7 i (k) 


estp (ik.R)J 

(k) “ cd^- i € 


(ilII-2) 


Here is sone function of k iiiiida does not involve R and .(k) 

My ~ ~ 3 

is a branch of the dispersion. The index 3 denoting the branch of the 
diuporsion is to be omitted for scalar models. The integration is cfver 
the fi ’ot Brillcfuin zone and £:->0 in the limit. For large R = j Rj 

tip integral can be expanded by the method of stationaiy phase -hhich has 

21 9,7,45,13 

bwOn discussed by Koster and subsequently used by many authors , 

Ti"'0 Gxpfinsion given below is a version of this method suitable for our 

purp'V’., . 

Uaing the formula 


X - it 


■bo 

i exp. £-i (x-i6 ) dt 


c ju.ition (AII- 2 ) can be written as 


^ ^ — I (R, ) 

-bo 


(AII-3) 


where 


I (R,'a) ) = ij ■iijjj exp^-i<jp(k,t)j (AlI-4) 

2 2 

with the phase (k,t)~^(0(k) - 


t - k.R 


(AII-5) 
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Wl^on Cp is stationaiy with respect to k aM t ™ has = 0 

, V7 <n ^ 

' kl ■" equations determine at least one (kj , t ) 

which specify the stationary point. Using (aII- 5) one finds that kj 

and t^ are such that 

o) . (k-^) = cO 
2 

and 


V 

k ' " 



E / t 
“ o 


(AII-6) 


Now one can expand (f (k,t) about this stationary point (k**,! ) 


O' o 


bj putting k = k" + K and t = t + ^ so that 

O Q 

£jp (k, t) = - kj.R + f^±^± t ^ ^i^ij^j (An-7) 


neglecting Iiigher order terms 






o 


and 


o 4 


ij 


dk.Tik. 


k = k 


(AII-8) 


(AlI-9) 


One c‘.an similarly expand >^(k) around k^ , bat only the^ (k^) term 
will /ti’vn a significant contribution because the higher terife will 
ratimatoly load to terns of the order of idiich will be negligible 

By 

wlion R is liirge. Further if there are more than one (kl, t^) that 
satisfy equation (AII-6) then a suinmatian must be done ewer all of them. 



108 .. 


Thus eqtie.tion (AII-4) can be inritten as 

1 (E, U ) = ZL e -0-- X (1,3) i ([( ^3^ exp-'C-r^fi fc 

^k3^t J JJJ 

— Or 0 ' 

+ 1 zi 


TMs can be further simplified by writing the phase in matrix and. vector 
notation as 


/V 


K. 

1 X 


+ t 



. T T A 

'tQK + t^KaiK 


\ifiiGre T" denotes the transpose. 

Further if we define p = K + -^ -i4_ 

o 



tlion 

-K + K oIk = -i ^ ■^*o£<^P 

0 ^ 

and d^ = d^p. We have used the fact that oC is symmetric which is 
obvious from its definition. Substituting this in equation (AH-IO), we 
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whore the limits of p integration is extended to ± Oo ^ Extending 
the limits of . integration similarly^ one obtains 







, . no 


• * 


Now for the scalar models, the branch index j can be omitted and 
comparing equation (AII-2) mth the corresponding equation of the scalar 
models described in equation (AI-1) one finds that X (k) is unity. 
Thus for the simple cubic, body centred cubic and face centred cubic 
scalar models one obtains 


C =- g (k^) 


3,B,F 



exp (ik^.R) 


all k 


R 


(AII-14) 


where evalixat mg Q _ and oC , as defined by equations (AII-8) and 
/ i y 

(AII-9) we got 


g (kj = 


S 4 Try 


2 2 2 
sin k + sin k „+ sin k 
■ 01 


q2 


o3 


2 2 -2 
L sin k cosk .cosk ^+sin k ^cosk ^cosk .+sin k cosk cosk 
01 o2 03 o2 o3 ol o3 ol o2 




(AII-15) 


£ (k^) = 


1 

8 TV 




(AII-16) 


where = cos ik . and S = sin . . 

1 01 1 Ol 


For face centred cubic 


g (k„) 


8irr 


2 2 2 2 2 2 


*Sx^s/s3^C^C2+C2C3+C3C,)Hti+C2)(C2+C3)(C3^-cp 


t 


2 2 2 
x-js^ {C2+C^)-^2 




Taii-17) 


For those scalar models, not more than one k is obtained as solution 

’ ~o 
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of equation (AII-6) in the limit of long waves because the constant 
frequency surfaces are splierical in this limit for these models. So 
the summation over is to be omitted. We have omitted this summation 
also in other cases where we are not interested in finding out numbers 
or exact res^fLts but require only qualitative aspects of the res-ults. 
For the case of the node! with polarization 

= eif(j) (k) e^^ (k) v/heTC (k) is the ^ th 

f* / / 

component of the jth polarization vector. Then one gets 




2 

id) ,R) 



(kj) (kj) exp(ikj.R) /R 

(AII-18) 


where 







The vector p emd the matrix OO in equation (AII-19) have to be 

evaluated at k^ . 

— o 



APPEWX III 


T^VALIT^TTON OF V - MATRICES 


Consider a lattice whose symmetry point group consists of elements 

iRf with the corresponding operations 0-r, . Let N be the order of 
C J « g 

„(^)J 

the group and let there be h irreducible representations f >V= 
of this group . tfe denote the dimensionality of the irreducible represent— 
ation j ' ' by n^; and the unitary matrices corresponding to operations 
in this irreducible representation by ^ ^ (R) whose traces or 
characters are given by X* '(H) . In addition to this point group 
of symmetry, the lattice has translational invariance but this is 
destroyed when we put in a substitutional impurity and ye assume that 
the system still has this point group symmetry. The displacenents of the 
substitutional impurity and its neighbours form a soace S Tdiich is 
invariant under the group of unitary operators • The characters 

^ (R) of this space S vrLth resnect to the group can be found and 
t^ien the representation of the group in the space S can be 
written as , 


4p 

where a 




n = 2ay p 

S R 


( A III -1 ) 
(A III ~ 2) 


This shows that the space S can be split into irreducible 
subspaces S ^ ^ ^ which occUr times for a given ^ . If there 

is a non-zero vector -p^in the space S then the operators, 


„ ( V) 


)*(K) 0, 


rs 


N 


g 


R 


rs 


R 


(A III-3) 


pro.iect out those vectors in the space S that belong to the r 


th 
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row of the V irreducible representation of the groupf^ This is 
subject to the condition that ^ is? not identically equal to 

03 ro . 


2 

. n 
>=1 


A( ) 

Since D ere n^ x n^^ matrices, there are in all 

^ = N such projection operators. Using these we can form the 

o 


different symmetiy adapted vectors out of . It is to be noted that 

} 

the choice of is completely arbitraiy and it is a matter of 
* i 

convenience. By vaiying r,s and ^ , we arrive at vectors that have 
the symmetries of different irreducible repiBsentations . These are 

orthogonal to one another and with proper normalisation they constitute 

A 

the columns of the recfuired Y matrix ' which block diagonolizes the 
matrices p and g . This block diagonolization follows the properties 
of these column vectors and from the fact that the operators 0„ 
comnute with p and g . 

’fe consider the specific case of the cubic group whose 

orcer N = 48 and whose character table is depicted in table 
g 

42 

III oelow in the notation of Hamermesh . Using equation (A III-3) 

for the simple cubic, body centred cubic , face centred cubic and 

the polarization model, which are denoted by the subscripts S, B, F 
and P respectively, one gets with the aid of table Til, 


n s = 2 hg " hu * 

r B = 2 ^Ig ^ J-lu * ^2g " V 

r'p-hu ® 

where @ denotes the direct product of the irreducible 
rep re sentpti ons . 
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TABLE III 

CHJSRACTBR TABLE FOR 



Hamennes 


:&or 


example ^ 


Wow following the procedure described we arrive 

Ax 


, ■ ' . at the folliarwing forms of V - mat rices. in these, for each 

irreducible representation of dimens-lonality appearing once in the 
reduced form of P , one has columns.- If it appears a^ 


times then one has a^ x n^ columns. 
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f \ 

The elements of IT-matrices given below are given in terms of 
numerical quantities given hy , 

b = (8)- ^ c = (24)" d = i , h = (12)“^, d' = (6)- i and 


h' = (2)' 




Vg = 


0 

0 

0 

0 

0 

V 

2A 


■Ig 


j cubic 

lattice 

we 

get 

/ 

0 

0 

O' 

0 

0 

0 

d' 

h' 

0 

0 

h 

d 

d' 

0 

h' 

0 

h 

-d 

d' 

0 

0 

h» 

-2h 

0 

d’' 

-h' 

0 

0 

h 

d 

d' 

0 

-h» 

0 

h 

-d 

df 

0 

0 -h' 

-2h 

0 

J 


lu 



E 

g 


(4 III-4) 


Here the row and column indices, starting from upper left comer are 
(0,0,0) , (1,0,0 ) , (0,1,0) , (0,0,1) , (-1,0,0) , (0,-l,0) and 
(0,0,-l) where the numbers refer to ( t ,m,n) values for the nearest 
neighbour sites according to equation (3.1) of chapter IH. 

(ii) Hor the body centred cubic lattice, we have 




(AIII-5 ) 


irith the row and column indices starting from theut>per left comer 
as (0,0,0), (0,1,0), (l,-l,0) , (0,-1, 1) , (-1,1,-!), (0,-l,0) , (-1,1,0), 
(0,l,-l) and (l,-l,l) , the numbers here referring to the ( t ,m,n) 
values for the nearest neighbour sites accord-ing to equations (4,1). 
(iii) and for the face centred cubic lattice. 
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with the row column indices starting from the upper left comer 

as (0,0,0) , (1,0,0) , (C,l,-1) , (0,1,0) , (l,0,-l) , (-1,1,0) ,(0,0,1) 

(-1,0,0) , (0,-l,l) , (0,-l,0), (-1,0,1) , (1,-1,0) and (0,0,-l) . 

A 

Finally the matrix Vp for the model with polarization is a 
27 X 27 matrix which can be partitioned into 9x9 submatrices as shown 
in equation (5.53) of chapter V. 


A 




(AIII-7) 


Fach of those submatitces have row and column indices the same as 
A 

doscilbed for Vg above ( equation (AIII-5) ). The older of columns 
A 

of Vp will coriespond to the following order of the irreducible 
representations of 0^ apoearing in the reduced form of . 


3F. , A, , E , 2F , A^ 

' Ig ' g ' 2g ' 2u 


lu 


, , F. and E 

’ 2u ' 1^ u 


This leads to the order of blocks described in chapter V. 


Explicitly, the matrices V^. ,(i,j = 1,2,3) of equation 
(ATTT-7) Can be written by using the folloviing nine dimensional 


vectors . 






T. = 
1 


= bd'- 


for i 


B}r writing these vectors side by side with aporopriate factors and 

A 

signs , the subisatricos V- . are, 

4- J 


V,, = (X, X„ X, 


X„ x^ X X ) 


LI '-1 "2 "0 0 0 8 0 0 7 ^ 

Vig = (^4 X^ X3 Xq X5 Xq Xq Xg ) 

^3 ^ ^^6 ^7 ^0 ^8 ^0 ^5 ^3 ^6 ^6 ^ 

''21 = (^0 ^0 ^ ^ ^2 ^0 ^0 ""O "^6 ) 

Vpp = (Y3 -X3 Z3 X^ Xq Xp Xg Xg Xq ) 

X -X X -X Z -X ) 

23 7 6 8 0 5 0 4 7 7 

'^31 = t^0 =^0 *7 ’'o ’'o "'e ^ ==2 ='0 ^ 

'^n K -X. x: -2Z„ X ) 

33 807634080 

■flien all these are put together in equation (ATII-7) one gets the 

A 

full 27 X 27 matrix Vp which is used in chapter V. 
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